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Abstract 

We continue the semiclassical analysis of the Loop Quantum Gravity (LQG) volume operator that was 
started in the companion paper |23| . In the first paper we prepared the technical tools, in particular the 
use of complexificr coherent states that use squares of flux operators as the complexifier. In this paper, 
the complexifier is chosen for the first time to involve squares of area operators. 

Both cases use coherent states that depend on a graph. However, the basic difference between the 
two choices of complexifier is that in the first case the set of surfaces involved is discrete, while, in 
the second it is continuous. This raises the important question of whether the second set of states has 
improved invariance properties with respect to relative orientation of the chosen graph in the set of 
surfaces on which the complexifier depends. In this paper, we examine this question in detail, including 
a semiclassical analysis. 

The main result is that we obtain the correct semiclassical properties of the volume operator for i) 
artificial rescaling of the coherent state label; and ii) particular orientations of the 4- and 6-valent graphs 
that have measure zero in the group S0(3). Since such requirements are not present when analysing dual 
cell complex states, we conclude that coherent states whose complexifiers are squares of area operators 
are not an appropriate tool with which to analyse the semiclassical properties of the volume operator. 
Moreover, if one intends to go further and sample over graphs in order to obtain embedding independence, 
then the area complexifier coherent states should be ruled out altogether as semiclassical states. 
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1 Introduction 



The volume operator in Loop Quantum Gravit}t^ (LQG) plays a central role in the quantum dynamics of 
the theory. Without it, it is not possible to define the Hamiltonian constraint operators [1], the Master 
constraint operators [5j, or the physical Hamiltonian operators [6j. The same applies to truncated models 
of LQG such as Loop Quantum Cosmology (LQC) [7J which is believed to describe well the homogeneous 
sector of the theory. In LQC the quantum evolution of operators that correspond to classically-singular 
Dirac observables typically remains finite. This can be traced back to the techniques used to quantize 
inverse powers of the local volume that enter into the expression for the triad, co-triad and other types of 
coupling between geometry and geometry, or geometry and matter [4j. Specifically, such derived operators 
arise as commutators between holonomy operators and powers of the volume operator. 

In view of its importance, it is of considerable interest to verify whether the classical limit of the LQG 
volume operator coincides with the classical volume. By this we mean that (i) the expectation value of the 
volume operator with respect to suitable semiclassical states which are peaked at a given point in phase 
space, coincides with the value of the classical volume at that phase space point, up to small corrections; 
and (ii) its fluctuations are small. 

It should be remarked that there are actually two versions of the volume operator that have been 
discussed in the literature [8l|9]. These come from inequivalent regularisations of the products of operator- 
valued distributions that appear at intermediate stages. However, only the operator in [9j survives the 
consistency test of [10], namely, writing the volume in terms of triads which then are quantised using the 
commutator mentioned above, gives the same operator up to /i-corrections as that obtained from direct 
quantisation. This consistency check is important as otherwise we could not trust the triad and co-triad 
quantisations that enter the quantum dynamics. 

A semiclassical analysis of the volume operator has not yet been carried out, although, in principle, 
suitable semiclassical (even coherent) states for LQG are available [11]. This is because the spectral decom- 
position (projection- valued measure) of the volume operator cannot be computed exactly in closed form. 
However, this is needed for exact, practical calculations. More precisely, the volume operator is the fourth 
root of a positive operator, Q, whose matrix elements can be computed in closed form |12j but which cannot 
be diagonalised analytically. 

In more detail, the volume operator has a discrete (that is, pure-point) spectrum and it attains a block- 
diagonal form where the blocks are labelled by the graphs and spin quantum numbers (labelling the edges of 
the graph) of spin-network functions (SNWF) [13]. SNWF form a convenient basis of the LQG Hilbert space 
[14j which carries a unitary representation of the spatial diffeomorphism group and the Poisson*— algebra of 
the elementary flux and holonomy variables [ISj. The blocks turn out to be finite-dimensional matrices whose 
matrix elements can be expressed in terms of polynomials of 6j symbols. Fortunately, these complicated 
quantities can be avoidecH by using a telescopic summation technique [16] related to the Elliot-Biedenharn 
identity [17], so that a manageable closed expression results. However, the size of these matrices grows 
exponentially with increasing spin quantum numbers and, since the expression for coherent states is a 
coherent superposition of SNWF's with arbitrarily large spin quantum numbers, a numerical computation 
of the expectation value using the numerical diagonalisation techniques, that are currently being developed 
[18] . is still a long way offH 

One way forward is to use the semiclassical perturbation theory developed in [19], and applied already 
in [201 121j . The basic idea is quite simple. In practical calculations one needs the expectation value of Q'^ 
where g is a rational number in the range ^ < q < \. In order to attain that, one has to introduce the 

^See dill] for recent books and [3] for recent reviews. 

^Notice that Racah's formula provides a closed expression for the 6j symbol but it also involves implicit sums and factorials 
that involve large integers. These quickly becomes unmanageable, even numerically. 

''The coherent superposition contains a damping factor that suppresses large spins, and thereby large matrices, so that the 
complicated infinite series over spin quantum numbers can be truncated at finite values, making only negligible errors. However, 
the computational effort required is currently too high, even for supercomputers; for example, see the estimates of computation 
time reported in [18] . 
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'perturbation operator', X := — 1, where the expectation value < Q > of the positive operator Q 
is exactly computable. Notice that X is bounded from below by —1. Then trivially < Q'^ >=< Q >^ < 
[1+X\^ >. Now we exploit the existence of positive numbers p such that the classical estimates 1+qx—px'^ < 
(1 + x)'' < 1 + qx are valid for all x > —1. Finer estimates of this form involving arbitrary powers of x 
are also available [20]. By virtue of the spectral theorem, this classical estimate survives at the quantum 
level and we have Y_ < Y < Y+ where Y+ = I + qX, Y_ = Y+ - pX^, Y = {l + X)i. It follows that 
<Y >e[<Y+> -p< X^ >, < Y+ >]. However, <Y+>=1 and < X^ >= [< > - < Q >^]/ < Q >^ 
is proportional to the relative fluctuations of Q, which are of order of h [llj. It follows that to zeroth-order 
in h we may replace < Q'' > by < Q which is computable, as are the error estimates in the way shown 
above. 

The exact computation of < Q >, < >, ..is feasible, but it is still quite involved. At this point 
it is convenient to use the fact that (see [H]), to zeroth-order in h, the computation of these expectation 
values (more generally, the expectation values of low-order polynomials in the flux operators) for SU{2) 
spin-network states coincides with the corresponding calculations for U{1)^ spin networks. On the latter, 
the volume operator is even diagonal. Hence we conclude that, as long as we are only interested in the 
zeroth-order in h contribution, we may evaluate the expectation value of the volume operator for a fictitious 
theory in which the non-Abelian group SU{2) is replaced by the Abelian group U{1)^. This dramatically 
simplifies all the calculations. 

The coherent states developed so far for LQG have all been constructed using the complexifier method 
reviewed in [22]. This generalises the coherent-state construction for phase spaces that are cotangent bundles 
over a compact group; see j24pl . This construction involves computing the heat-kernel evolution of the delta 
distribution (which is the matrix element of the unit operator in the Schrodinger (position) representation) 
with respect to a generalized Laplace operator, called the 'complexifier', followed by a certain analytic 
continuation. 

Now the unit operator in LQG can be written as a resolution of unity in terms of SNWF's and, although 
the heat kernel is a damping operator, since the SNWF are not countable (the LQG Hilbert space is not 
separable), the resulting expression is not normalisable. However, it does give a well-defined distribution 
(in the algebraic dual of the finite linear span of SNWF's) which can be conveniently written as a sum 
over cut-off states labelled by finite graphs. These states, called "shadows" in [24], are normalisable, and 
can be used to perform semiclassical calculations. Of course, one expects that the good semiclassical states 
will only be those cut-off states that are labelled by graphs which are sufficiently fine with respect to the 
classical three-metric to be approximated. 

Thus the input in the semiclassical calculations consists in the choice of a complexifier and the choice of 
a graph. One may wonder why the complexifier has to be chosen a priori rather than being dictated by the 
dynamics of the theory: that is, why do the coherent states remain coherent under quantum time evolution? 
For example, for the harmonic oscillator, the complexifier is just the Laplacian on the real line. 

The reason is two-fold. On the one hand, one could perform a constraint quantisation so that we 
are working at the level of the kinematical Hilbert space on which the quantum constraints have not yet 
been imposed. In this situation, all that is needed is to ensure that the volume operator, which is not 
gauge invariant (and thus does not preserve the physical Hilbert space), has a good classical limit on the 
kinematical Hilbert space. We would certainly not trust a constraint quantization for which even that was 
not true. 

On the other hand, one could also work at the level of the physical Hilbert space, as outlined in [6]. 
Then, one would expect that the time evolution of any reasonable choice of complexifier coherent states 
with respect to the physical Hamiltonian should remain coherent (and peaked on the classical trajectory) 
for sufficiently short time intervals. Note that for interacting theories, even the simple example of the 
anharmonic oscillator, globally stable coherent states have yet to be found. 

The choice of the complexifier will be guided by practical considerations: namely (i) it is diagonal on 
SNWF's; and (ii) it is a damping operator that renders normalisable the heat-kernel evolution of the delta 

''See also [25| for related ideas valid for Abelian gauge theories such as Maxwell theory or linearised gravity. 
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distribution restricted to a graph. Moreover, it should be gauge invariant under the SU{2) gauge group. As 
for the choice of graph, for practical reasons one will use graphs that are topologically regular, that is, have 
a constant valence for each vertex. Indeed, the semiclassical calculations performed in [20\ [2T] were done 
using graphs of cubic topology, with good results. 

The existing literature on such coherent states for LQG can be divided into two classes, depending on 
a certain structure that defines them. The first are gauge-covariant flux coherent states, which depend 
on collections of surfaces and path systems inside them [11]. The second are area coherent states, which 
also depend on collections of surfaces but involve area operators rather than flux operators [22]. In our 
companion paper [23j we have reviewed and generalised both constructions. In its most studied incarnations, 
the collection of surfaces involved in pjj is defined by a polyhedronal partition of the spatial manifold, while 
the collection used in [22] uses a plaquette of foliations of the spatial manifold. 

The main objective of this series of two papers is to analyse the semiclassical properties of the volume 
operator with respect to both classes of coherent state. In [23j we employed the states of [llj. In this paper 
we employ the states of [22] . This exhausts all known coherent states for LQG. The result of our analysis 
is that, if we use the states [11], the correct semiclassical limit is attained with these states for n = 6 only. 
If instead we use the states [22], the correct semiclassical limit is attained only for: 1) artificial rescaling 
of the coherent state label; and 2) particular embeddings of the 4-valent and 6-valent graphs with respect 
to the set of surfaces on which the complexifier depends. However, the combinations of Euler angles for 
which such embeddings are attained have measure zero in SO (3), and are therefore negligible. Thus the 
area complexifier coherent states are not the correct tools with which to analyse the semiclassical properties 
of the volume operator. 

If one wants to obtain embedding independence, a possible strategy is to sample over graphs (Dirichlet- 
Voronoi sampling [H]) as outlined in [11]. What this strategy amounts to is that, instead of singling out 
one particular coherent state ip^y^m — as defined in terms of a single graph 7 — one considers an ensemble of 
coherent states constructed by averaging the one-dimensional projections -Py,m onto the states ^^7,™ over a 
subset Tm of the set of all allowed graphs. In other words, one considers a mixed state (with an associated 
density matrix) rather than a single coherent state. In such a way, if the subset Tm is big enough, it can be 
shown ([11]) that it is possible to eliminate the embedding dependence (the 'staircase problem^). 

It is straightforward to deduce that the area complexifier coherent states cannot be used to construct 
embedding-independent, mixed coherent states because of condition 2) above. We thus claim that area 
complexifier coherent states should he ruled out as semiclassical states altogether if one wants to attain 
embedding independence. Instead one should use the fiux coherent states, as was done in [llj. For such 
states we will show that the correct semiclassical limit is attained only for n = 6 . In other words, up to 
now, there are no semiclassical states known other than those with cubic-graph topology! 

Thus the implication of our result for LQG is that the semiclassical sector of the theory is spanned 
by SNWF that are based on cubic graphs. This has some bearing for spin- foam models |27j which are 
supposed to be — but, so far, have not been proved to be — the path-integral formulation of LQG. Spin 
foams are certain state-sum models that are based on simplicial triangulations of four manifolds whose dual 
graphs are therefore 5-valent. The intersection of this graph with a boundary three-manifold is 4-valent, and 
therefore we see that spin-foam models, based on simplicial triangulations, correspond to boundary Hilbert 
spaces spanned by spin- network states based on 4-valent graphs onljj^]. However, we have proved that the 

^Roughly the staircase problem can be stated as follows. Consider an area operator As for a surface S. If we compute the 
expectation value for As with respect to a coherent state ip-y,m, such that the surface 5* intersects transversely one and only one 
edge e of 7, then the expectation value of the area operator coincides with the classical value A{m). However, if the surface S 
lies transversally to the edges, then we do not obtain the correct classical limit. 

^ As an aside, whether this boundary Hilbert space of spin foams really can be interpreted as the 4-valent sector of LQG is 
a subject of current debate, even with the recent improvements [28) in the Barrett-Crane model [29]. There are two problems: 
first, the boundary connection predicted by spin foams does not coincide with the LQG connection [3D]. Secondly, the 4- 
valent sector of the LQG Hilbert space is not a superselection sector for the holonomy flux algebra of LQG. In fact, the LQG 
representation is known not only to be cyclic but even irreducible [31]. Therefore the 4-valent sector is not invariant under the 
LQG algebra. 
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correct semiclassical states for analysing the semiclassical properties of the volume operator are the gauge 
covariant flux states. For such states, only those of cubic topology give the correct semiclassical value of the 
volume operator. 

Even if the mismatches between the 4-valent sector of LQG and the boundary Hilbert space of spin 
foams could be surmounted, the result of our analysis seems to be that the boundary Hilbert space of current 
spin-foam models does not contain any semiclassical states! This apparently contradicts recent findings that 
the graviton propagator derived from spin-foam models is correct [32j. However, it is notable that these 
latter results only show that the propagator has the correct fall-off behaviour: the correct tensorial structure 
has not yet been verified. 

One straightforward way of possibly repairing this situation is to generalise spin-foam models to allow 
for arbitrary — in particular, cubic — triangulations, as suggested in [331 |3l]. 

In the present paper we will perform detailed calculations of the expectation value of the volume operator 
for the 4-, 6-, and 8-valent graphs using the area coherent states reviewed in |23) . The main difference 
between the flux coherent states and the area coherent states is that the former depend on a discrete set of 
flux operators while the latter depend on a continuous family of area operators. In particular, the surfaces 
on which the area operators depend fill all of space. Hence the question arises whether the area coherent 
states depend less severely on the relative orientation of the chosen cut-off graph with respect to the system 
of surfaces. Specifically, it would be interesting to know if the calculations of the expectation values are 
insensitive to Euclidean motions of the graph (translations and rotations) when the spatial metric to be 
approximated is the Euclidean one. 

The structure of the paper is as follows. In Section [2] we calculate the expectation value for a general 
graph. The expectation value depends on the edge metric, already encountered in [23], which is close to 
diagonal when the edge lengths, 5, are much larger than the lengths, I, of the boundaries of the surfaces. 
In order to compute the off-diagonal metric entries we stick to the graphs defined in our companion paper 
as otherwise we encounter severe 'bookkeeping' problems. These graphs are dual to tetrahedronal, cubical 
and octahedronal triangulations, all of which derive from a cubulation of the spatial manifold, as discussed 
in [23]. One can expand the expectation value into a convergent power series in 1/6. We establish that, 
contrary to the results in [23], to zeroth-order in both h and 1/5 the correct expectation value is obtained 
for n = 6 and n = 4, but only for specific embeddings of the graph with respect to the surfaces on which the 
area operators depend. However, we show that such embeddings have measure zero in SO{3). Moreover, 
by using the coherent states in [11] we also find that the graphs (4-, 6, and 8-valent) are all Euclidean 
invariant, up to sets of measure zero in S0{3) and up to linear order in 1/5. Nonetheless, the embeddings 
with non-zero measure all fail to give the correct semiclassical value for the volume operator. 

This shows that coherent states |llj are not the correct tools for analysing the semiclassical properties 
of the volume operators. 

In Section 3 we perform an explicit calculation of the expectation value of the volume operator for 
4-, 6-, and 8-valent graphs, including the higher-order effects in 1/5. This contribution is only Euclidean- 
invariant up to linear order. For each of these valences we first consider a graph constructed in terms of 
regular simplicial, cubical and octahedron cell complexes, as done in [23]. After that, we study the effects 
of rotations and translations. 

In Section 4 we summarise our results and draw appropriate conclusions. 

The more technical calculations have been transferred to an appendix in order not to distract from the 
line of argument in the main text. 

2 Volume Operator Expectation Values for Area Coherent States 

In this Section we compute the expectation value of the operator V^^„ for an arbitrary n-valent vertex, v, 
for the stack family coherent states using the replacement of SU{2) by U{1)^. This uses the calculational 
tools developed in [23J, to which we refer for our notation. We may, therefore, replace the SU{2) right- 
invariant vector fields by U{1)^ right-invariant vector fields ,s = ih'jd/dh-^j acting on /i]- := A^{ei{v)). 
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The crucial simplification is that these vector fields mutually commute. Their common eigenfunctions are 
the spin-network functions which, for U{1)'^, take the explicit form 



eeE(7) j=l 



(2.1) 



We will refer to them as 'charge network states' because the G Z are integer valued. Using the spectral 
theorem we may immediately write down the eigenvalues of Vy^^ on T.y^„ as 



\ — 



1 



nf' rf/ nf^ 



l</< J<i<'<JV 



(2.2) 



What follows is subdivided into four parts. We begin by performing the calculation for a general graph. 
This leads to the inverse of the edge metric which, for large graphs, is beyond analytical control. In the 
second part we restrict the class of graphs, which lets us perform perturbative computations of the inverse 
of the edge metric. This gives a good approximation of the actual expectation value. In the third and fourth 
parts we consider the dependence of our results on the relative orientation of the graph with respect to the 
family of stacks. 



2.1 General expectation value of the volume operator 

In this Section we drop the graph label and set t'^^, := 5^^ ^Ce" t := Pp/a? . This gives a positive, symmetric 
bilinear form on vectors n := {'nfj)e&E{-y), j=i,2,?, that is defined by t{n,n') := X^e e' j fc ^j'^'fc '*ee' = ■ t ■ n' . 
We also set ■ n := Yle j Z^{e)nj. 

The coherent state associated with an n-valent graph in which more than one edge intersects a given 
plaquette S, is as follows: 

^^,7 = E e-^""*-" e^"-r^,„ (2.3) 



The norm of the coherent states is given by 



E 



(2.4) 



where 



P^{e)=i j{Z-A) =: ^E] 



(2.5) 



The length parameter, b, that appears here is generally different from the parameter, a, that enters the 
classicality parameter t = 1"^, / a? , as explained in [23] . 
The expectation value for the volume operator is 



<V„> 



V >Z,7 = 



'2,7 1 



2,7 1 



where 



A7,.(n)=4./|^ E ^(^'^''^") det 



n 



(2.6) 



(2.7) 



er\e'r\e"=v 

are the eigenvalues of the volume operator V^^^. We have introduced the notation 



det (n) := 



7 



The semiclassical limit of the volume operator is obtained from (j2.6|) in the limit of vanishing t. That 
is, it is the zeroth-order in t of the expansion of (j2.6p in powers of t. Since (j2.6p converges slowly for small 
values of t, we will perform a Poisson transform which replaces thyj which converges quickly. To this end, 
analogous to [20], we introduce the following variables 



tee J 



Te := VTe, x) := T,n% y][x) := x^T,, C| := Pj/Te, w^{n) := n'^j/T^, Ai^, := A,,,6^' := 

(2.9) 

Notice that the diagonal entries of A equal unity. The off-diagonal ones, however, are bounded from above 
by unity by the Schwarz inequality applied to the scalar product defined by A and are restricted to the 
six-dimensional subspace restricted to vectors with non-zero entries for the e, e' components only. 
Then (j2.6p turns into 

E„ e--^-^-^ e^'^-"A-, ,,(j/(x)) 
Applying the Poisson transform to (j2.10p we obtain (recall N := |£'(7)|) 



< >^'7- V- .-xT.A.. oc-x (2-10) 



. _ E„ /m3^ .^^^xe-^— (")^-e-^-^--e^^-A,,.(y(x)) 

< >Za ^ j^^^ ^3jy^ ^-2mw{nY-Xf^-xT^-A-x ^2C-x \^-^^) 

In order to perform the Gaussian integrals in (j2.1ip we notice that, by construction, ^ is a positive- 
definite, finite-dimensional matrix, so that its square root, and its inverse are well defined via the 
spectral theorem. Hence we introduce as new integration variables 

z| := Y.^VA)f;xi, y]{z) = ^x] = ^Jl {\^r%4 (2.12) 

e\k e'k 

The Jacobian from the change of variables drops out in the fraction, and ()2.1ip becomes 



Now one would like to shift z into the complex domain by v yl {C — it^w) and then perform the ensuing 
Gaussian integral. This is unproblematic for the denominator of (12. lip which is analytic in z; however, the 
numerator is not. The careful analysis in [20] shows the existence of branch cuts in of the fourth-root 
function involved. In turn, this shows that, in the semiclassical limit, both numerator and denominator are 
dominated by the n = term while the remaining terms in the series are of order h°° (i.e they decay as 
exp(— /c„/t), t = £p/a'^ for some kn > 0, lim„^oo kn = oo). See [20] for the details. 

The upshot is that to any polynomial order in h we may replace (I2.13P by 

" ^'^^ r — d^^z e-^- ^ ^ ^ 

J]R3iv a z e 

which is now defined unambiguously because the argument of the fourth root is the square of a real number. 

The denominator of ()2.14p simply equals -y/vr . Therefore, the only ^-dependence of (j2.14p lies in the 
numerator in the function A^^^,. We now note that the eigenvalues of the volume operator come with a 
factor of ^p, as displayed in (|2.7p . Pulling it under the square root into the modulus, and noticing that the 
modulus is a third-order polynomial in y, we see that 

A^,,(y(z + VI^'C)) = . /|^ E i^py(^ + -^''C)) I (2.15) 

V ene'ne"=v 
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where we have 



e',k 
e',fc 

where ()2.5p has been used in the last hne. 

To extract the leading order in t of (j2.14p is now easy. First note that the matrix elements of both A and 
Tf./\/t are of order unity. Since a is some macroscopic length scale, the first term that is proportional to z in 
(j2.16p is therefore of order -v/t, while the second is of zeroth-order in t. Therefore F := X^^y(y{z + \/]4 C)) 
is of the form 



F{zVt) = {lQ + P{zVt) (2.17) 

where P is a certain sixth-order polynomial in z^/t with no zeroth-order term, while Q is independent of z. 
Moreover, Q + P{zVt) is non-negative for all z because it is the square of a third-order polynomial in z. In 
particular, this holds at z = 0, and therefore Q is also a non-negative number. Then, provided Q > 0, we 
can define 

f{zVt) := =: ^l + R{zVt) (2.18) 

where i? is a sixth-order polynomial with no zeroth-order term which is bounded from below by —1. Now, 
as in [19] , we exploit the existence of r > such that 

1 + -R-rR^ < f < 1 + -R 2.19 
4 - ■' - 4 

for all i? > —1. Inserting this estimate into (I2.14P we can bound the integral from above and below because 
the Gaussian is positive. The integrals over R and R^ are finite and are at least of order t because odd 
powers of z do not contribute to the Gaussian integral. 
It follows that to zeroth-order in t we have 



48 

ene'ne"=v 



<K>z,,= j^ E e{e,e',e") ^detjY) (2.20) 

where 



v--=(lYE^i^-']r^E'^ (2-21) 



bJ ^-^ TpTp' 

e'k 



This is as far as we can go with the calculation for a general graph. Notice that the inverse of the edge 
metric appears in this expression and, for a general graph, this is beyond analytical control. Therefore we 
will now make restrictions on the graph so as to analyse (|2.2ip further. 



2.2 Simplifying assumptions 

The assumptions about the class of graphs to be considered are as follows: 
1. Coordinate Chart 

The graph, the region R and the families of stacks lie in a common coordinate chart X : — > a. This 
is not a serious restriction because the general situation may be reduced to this one by appropriately 
restricting attention to the various charts of an atlas that covers a. 
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2. Tame Graphs 

We assume that the graph is tame with respect to the stacks. By this we mean that for each direction 
/, and each stack a, a given edge, e, of the graph enters and leaves that stack at most once. This 
means that the graph does not 'wiggle' too much on the scale of the plaquettes. Analytically, it means 
that /^"^ vanishes whenever lA'^e] > 2 for any e, and that, for given e, the number l^j^^ is non- vanishing 
at most for either A'^e = +1 or Ng = —1 but, not both, and independently of a. Finally, it means that 
the sets S^j^'^ are connected. 

3. Coarse Graphs 

We assume that the graph is much coarser than the plaquettation in the sense that any edge intersects 
many different stacks in at least one direction /. 

4. Non-Aligned Graphs 

We exclude the possibility that distinct edges are 'too aligned' with each other in the sense that the 
number of stacks that they commonly traverse is much smaller than the number of stacks that they 
individually traverse . 

Pictorially, the situation therefore typically looks as in figure [TJ A consequence of the tameness, coarseness 




Figure 1: Example of a tame, course and non-aligned graph 

and 'alignedness' assumption is that \tl^,\ <^ tj = tJe, tj, = tj,^, for all e,e', as is immediately obvious 
from the formulae displayed in our companion paper because the number of stacks with |A^e| = |A^e/| = 1 
will be very much smaller than the number of stacks with |A^e| = Ij = or |A^e'| = Ag = 0. Hence 
the edge metric will be almost diagonal. This is important because we need its inverse, which can only 
be calculated with good approximation (that is, for large, semiclassically relevant graphs) if it is almost 
diagonal. The graphs that we will eventually consider are embeddings of subgraphs dual to tetrahedronal, 
cubical or octahedronal triangulations of M^. These correspond to embeddings of regular 4-,6-,8-valent 
lattices, which ensures the non-alignedness property. 

Thus, without loss of generality, we may choose the stacks and plaquettes as follows. Using the availability 
of the chart X : — > fi; s ^ ^(s) we consider the foliations defined by the leaves Ljt := X/(M^) 
where for ejjx = 1 we set Xl{v},v?) := X{s^ := t, s'^ := u^,s^ := u^). The stacks are labelled by 
a = (a^,a^) G Z-^. The corresponding plaquettes are given by pi^ = {A/([q + u]l); u £ [0, l)^} where / > 
is a positive number. 

Likewise, using the availability of the chart, we take the edges of the graph to be embeddings of straight 
lines in (with respect to the Euclidean background metric available there), that is, e{t) = X{se + Ve6t) 
where Ve is a vector in and e(0) = X{se) defines the beginning point of the edge. 
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After these preparations, we can now analyse (j2.2Up and (|2.21|) further. Recall that 

^i=E/ dtE,{pf)a{pf,e) (2.22) 

and 

*ee' = E / dta{pf,e)a{pf,e') (2.23) 

al 

By the assumption about the graphs made above, the signed intersection number takes at most the numbers 
±1 and independently of a, so that cr(pf^,e)^ = ala{pf^,e) for certain = ±1 which takes the value +1 
if the orientation of e agrees with that of the leaves of the foliation, —1 if it disagrees, and if it lies inside 
a leaf. If we assume that the electric field is slowly varying at the scale of the graph (and hence at the 
scale of the plaquettes as well) then we may write 

E]^Y.*2AeM) (2.24) 

where pi = p'^^l^^-j^ and v is the vertex at which e is adjacent and which is under consideration in Vy^y. It 
follows that ()2.2ip can be written as 



e' ,k I e' ' I 

where we have used = ttJe- 

Now, by construction, A = 1 + B with B off'-diagonal and with small entries 

Bee' = -r== (2-26) 

which are of the order of 1/6 since two distinct edges will typically only remain in the same stack for a 
parameter length /, while the parameter length of an edge is 6. Now notice that under the assumptions we 
have made, we have Ij^, = if e, e' are not adjacent. Define Se to be the subset of edges which are adjacent 
to e. Then 

11-^^11^ = ^ ^ Xe'Be'eBee"Xe" 
e e',e"eSe 

l2 



< [snpBl,] ^ [ ^ xe'^ 

£ (^)'E [(Ei=)"'(E4)"^ 
£ (s)'*'E E- 



6Se e'eSe 

\ )\/r \ \ ™2 



e' e 

< (:^y^2||^||2 (2.27) 



.5. 

Here, in the second step we estimated the matrix elements of B from above; in the third step we applied the 
Schwarz inequality; in the fourth step we estimated jS'el < M, where M is the maximal valence of a vertex 
in 7; and in the sixth step we exploited the symmetry 

XsAe) = \ ^ ^,^^^0 =Xsde) (2.28) 
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as well as the definition of the norm of x. 

It follows that for 1/5 < M, B is bounded from above by unity. Therefore, the geometric series = 
1 + ^)"' converges in norm. Hence we are able to consider the effects of a non-diagonal edge metric 

up to arbitrary order, n, in 1/5. Here we will consider n = 1 only and write A^^ = 1 — {A — 1) = 2 ■ 1 — A. 
But before considering corrections from the off-diagonal nature of A notice that, to zeroth-order in 1/5, 
equation (j2.25p becomes simply 

aiE.ipi) (2.29) 



ye 

3 



Inserting (p:29]) into ([2:20|) we find 



dei{E){v)\ J ^ <e,e',e") ^det^/a) det(p, 



ene'ne"=ii 



e,e ,e" 



(2.30) 



where 



det(p,) = ^.euKe"'"' 



[0,1)2 



[0,1)2 



(2.31) 



On recalling that </ = eabcX^^u^X^c^^^^^ with X^?(ii) = X'^(s^ = = + /u^ = + /'u2) for 
e/jii" = 1, we find 



Hence (|2.30p becomes 



det(p^) « f[det{dX{s)/ds)\ 



X{s)=v 



(2.32) 



<K>z,7~ X Vl^^t(^)(i;)| [det(5X(s)/as)] 



I X(s)=t) 




e(e,e',e") 



det 

e,e',e" 



(2.33) 



We can draw an important conclusion from expression (j2.33p . Namely, the first three factors approximate 
the classical volume Vv{E) as determined by E of an embedded cube with parameter volume {al/b)^ . When 
we sum (|2.33p over the vertices of 7, which have a parameter distance, 5, from each other where / <C 5 by 
assumption, then the volume expectation value only has a chance to approximate the classical volume when 
the graph is such that 5 = al/h, or 5/1 = a/h. This could never have been achieved for h = a and explains 
why we had to rescale the labels of the coherent states by (a/b)^ while keeping the classicality parameter at 
t = tp/a^. See our companion paper for a detailed discussion. There we also explained why one must have 
5/1 actually equal to a/6 and not just of the same order: While one could use this in order to favour other 
valences of the volume operator, the expectation value of other geometrical operators such as area and fiux 
would be incorrect. 

Assuming 5/1 = a/h we write (I2.33P as 



< K > 



V,iE) G 



(2.34) 



thereby introducing the graph geometry factor G^^^. It does not carry any information about the phase 
space, only about the embedding of the graph relative to the leaves of the three- foliations. From the fact 
that (I2.34P reproduces the volume of a cube up to a factor, we may already anticipate that the geometry 
factor will be close to unity for, at most, a cubic graph. Whether this holds for an arbitrary orientation of 
the graph with respect to the stack family will occupy a large part of the analysis which follows. 



2.3 Preliminary analysis of the graph geometry factor 

We start by investigating the behaviour of the graph geometry factor G^^y under diffeomorphisms, (p, of a, 



that is, under G. 



7,1; 



while the linearly-independent families of stacks are left untouched. This 
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will answer the question of how much the geometry factor depends on the relative orientation of the graph 
with respect to the stacks. 

In fact, the orientation factor e(e,e',e") is invariant under diffeomorphisms of the spatial manifold a. 
The signature factor 



is obviously invariant under any diffeomorphism that preserves the foliations F , i.e. which map leaves onto 
leaves, because 



where Lu is any leaf in t which intersects e transversely. Since we consider graphs whose edges are embedded 
lines in with the same embedding that defines the stacks, it follows that the geometry factor is invariant 
under any embedded global translations in R^. 

Next, since global rescaling in M'^ preserves the foliations and the topological invariant (j2.36p . the geom- 
etry factor is also invariant under embedded global rescalings of M"^. Finally, any embedded global rotations 
of M'^ that preserves all the orientation factors a], will leave the geometry factors invariant. Since the ori- 
entation factors only take the values +1,-1,0 (depending on whether an edge agrees, disagrees with the 
orientation of the leaves, or lies within a leaf), there will be a vast range of Euler angles for which this con- 
dition is satisfied, if the graph is an embedded, regular lattice of constant valenc^. Hence, in order to check 
whether the geometry factor is rotationally invariant under any rotation we need only worry about those 
rotations which lead to changes in the a],. Likewise, if we rotate a graph which is dual to a polyhedronal 
complex, we expect that the expectation value remains invariant as long as the graph remains dual to the 
complex. 

Fortunately, using the explicit formulae derived for the edges and vertices for n = 4-, 6-, 8-valent graphs 
displayed in [23] we can calculate the o], for each edge e. Intuitively, it is clear, that whenever many of the 
(Tg change from +1 to —1, we can expect a drastic change of the expectation value. However, one has to take 
into account the combined effect of these changes, and this is what makes rotational invariance possible. As 
a first step we determine the action of a rotation on the sign factors. 

2.4 Explicit calculation of the erf terms 

In what follows we discuss the cases that show a drastic change in the value of det^ ^ (cr) = e/ji^-CTgCJ^/Cr^,, 
caused by a change in the values of CTg. To carry out this calculation we will perform a rotation of each of 
the three different types of lattice analysed so far: namely, the 4-, 6- and 8-valent lattices. These rotations 
will be parameterised by Euler angles and will be centred at a particular vertex of the lattice, for example 
Vq. The effects of a rotation will depend on the distance of the vertices from the centre of the rotation. 
In fact, the position of each vertex in the lattice after rotation will depend on both the distance from the 
centre of the rotation and the Euler angles used in the rotation. Fortunately, the values of the terms a\ will 
not depend on the former but only on the latter. 

This is easy to see since the value of fig can be either 1, -1 or depending on whether the edge is 
outgoing, ingoing, or lies on the plaquette in the direction I. Thus it will only depend on the angle the edge 
makes with the perpendicular to the plaquette in any given direction, i.e. it will depend on the angles the 
edge makes with respect to a coordinate system centred at the vertex at which the edge is incident. Clearly, 
only the values of the Euler angles of the rotation will affect the angles each edge has with respect to the 
vertex at which it is incident. In particular, since the graph we are using is regular, following the rotation, 
all edges which were parallel to each other will remain such, and thus will have the same angles with respect 
to the vertex at which they are incident. This implies that in order to compute the values of the terms a^, 
we can consider each vertex separately and apply the same rotation to each vertex individually. 

^In fact, for a random graph we may also have rotational invariance on large scales. 




(2.35) 




(2.36) 
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On the other hand, the distance from the centre of the rotation affects the position of each vertex with 
respect to the plaquette structure, and thereby affects both the values of the terms t^^i and the number of 
them that are different from zero. These effects can be easily understood with the aid of the two-dimensional 
diagram (Figure [2]). 



e 




Figure 2: Example of translation and rotation 



It is clear that, for any two parallel edges, the angle each of them has with respect to the vertex at which 
they are incident, is independent of the distance of the edge from the centre of rotation. On the other hand, 
the values of the t^^i will depend on both the rotation and the distance of the centre of rotation, since the 
position of the rotated vertex with respect to the plaquette depends on both these parameters. Therefore, 
we can tentatively assume that two different geometric factors will be involved in the computation of the 
volume operator: 

i) G-yy, which indicates how the terms are affected by rotation. This geometric factor affects all 
orders of approximation of the expectation value of the volume operator. 

ii) C^y, which indicates the effect of rotation on the terms t^^i. This term affects only the first- and 
higher-order approximations of the expectation value of the volume operator, not the zeroth-order. 

In what follows we will analyse the geometric term G^y: i.e., we will analyse the changes in the values 
of the (Tg due to a rotation applied at each vertex independently. We will do this for the 4-, 6- and 8-valent 
graphs separately. The geometric factor C^y will be analysed in subsequent Sections. 

As we will not see, our calculations show that for all 4-, 6- and 8-valent graphs, the rotations that 
produce drastic change in the values of det^ (a) = euKcr^o-^iaK have measure zero in 50(3) since they 
occur for specific Euler angles rather than for a range of them. 

Let us start with the 6-valent graph (Figure [3]). From the discussion above, we need only consider the 
effects of the rotation on one vertex, Vq. 

In order to compute the change in the values of the individual Ug, we will divide into eight small sub- 
cubes the cube formed by the intersection of the plaquettes in the three directions and containing the vertex 
we are analysing (Vq). It is then easy to see that for each edge e, the corresponding value of depends on 
the sub-cube in which it lies. In particular, we have the following table for the values of a^. 
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Figure 4: Division of the cube in 8 sub-cubes 
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From the above table it is clear that when an edge moves from one of the eight cubes to another, the 
values of each of the changes accordingly. Given any 6-valent vertex, each of the six edges incident at 
a vertex will be in one distinct cube. Moreover, since any two edges incident at a vertex can be either 
co-planar or perpendicular (in the abstract pull-back space with Euclidean metric), there are only certain 
combinations of allowed positions. For instance, for the edges ei, 62, 63, 64, 65, eg only the combinations 

CiB2DsEiF5HQ , CiD2AsFiG5Ee , AiB2CiHiE^GG , BiA2DiHiG^FQ (2.37) 

are allowed (here the notation Ai means that the edge i lies in the cube A); the combination AiB2C^DiG^FQ 
is not allowed. 

Because of the highly symmetric structure of the 6-valent graph we do not have to analyse all possible 
combinations of all the six edges incident at a vertex, since different combinations are related by symmetry 
arguments. For example, the combination in which edges ei, 62, 63, 64, 65, eg lie in the cubes GiD2A^F4,G^Eq^ 
and the combination in which they lie in the cubes FiG2H^GiiD^AQ, lead to the same value of | det^ ^1 ^1 ((t)| = 
k/jii'<7p0"'^,cj^, I, and an equal number of det / „"(ct) > and det„ ' „"(o") < 0, but obtained from different 
triplets e, e , e . In particular, any consistent relabelling of the edges will produce the same overall result 
for the determinants of the triplets. These symmetries reduce considerably the number of cases that need 
to be analysed. 

In what follows, we consider the cases for which the edges ei, 62, 63, 64, 65, eg lie in the following combi- 
nations of cubes: 

GiB2D^E^F^Hq , GiD2AsFiG5Ee , AiB2GiHiE^Gf, , BiA2DiH^G5F^ (2.38) 

For each of these cases there will be sub-cases according to whether one edge or more lie in a particular 
plaquette, or are parallel to a given direction /, J, K. These sub-cases are the following: 

1. No edge lies in any plaquette, or is parallel to any of the directions. 

In this case we obtain | det^ g/ g" (o")| = 4 for all triplets, but four of these triplets will have det^ ^1 ^1 (a) = 
—4 while the remaining four will have det / // (a) = 4. 

^ e,e ,e V 

2. Only one edge lies in a particular plaquette (say the J direction) (see Figure [3l). This edge and its 
co-linear edge will have equal to zero (J being the direction of the plaquette in which the edge lies.) 

In this case we obtain det^ ^' ^" (a) = —4 for four triplets, and det^ ^1 ^" (a) = 4 for the remaining four 
tripletsl!| 

3. Two edges lie in two different plaquettes such that each of these two edges and their respective co- 
linear edges will have equal to zero in the direction of the plaquette in which they lie. In this case, 
because of the geometry of the 6-valent lattice, the remaining edges will each be parallel to a given 
direction J such that all but the are zero. 

In this case we obtain det^ ^1 ^1 (cr) = 2 for four triplets while the remaining four will have det^ ^ ^1 (o") = 
-2. (See Figure ED. 

4. Each edge is parallel to a given direction such that all the a], (for any /, J, K) are equal to zero except 
for the one in the direction to which the edge is parallel. In this case we obtain detg ^ ^1 (fi) = 1 for 
four triplets and det^ ^ e"('^) = — 1 for the remaining four. (See Figure [6]) 

*Note that the geometric factor associated to this edge orientation will coincide with the geometric factor as derived from 
case 1). In this sense, case 2) can be seen as a limiting case of 1) 
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Only sub-cases 3 and 4 might lead to a change of value for the geometric factor G^y. However cases 2, 3 
and 4 have measure zero in 5*0(3). 

As a demonstrative calculation on how this is derived we will choose case 3. In particular, we select the 
configuration depicted in Figure [5] which can be obtained by a rotation of the original configuration in Figure 
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Figure 5: Example of configuration with zero measure in 50(3) 
El Let us consider the linearly-independent triples comprised of the edges that connect the barycenter of 
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Figure 6: Regular 6-valent graph 



the cube to the vertices Vi, V2 and V3. In the original configuration, the coordinates of these vertices are 
(in what follows we will denote the length of an edge, e, by (5e = (5) 

Vx = (0, 0, 6) ,V2 = (5, 0, 0) , ^3 = (0, <5, 0) (2.39) 

By applying a general Euler rotation, whose matrix representation is given in (j3.2.2p . the coordinates of the 
rotated vertices become: 

Vi = {R135, R23S, R335) , V2 = {Rii5, R21S, R31S) , V3 = {R12S, R226, R32S) (2.40) 

Our task now is to determine which Euler angles would give rise to the configuration in Figure [5j Since 
in such a configuration the edges eoi (the edge joining the barycenter of the cube to vertex Vi) and eo3 lie 
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in the plane x-y, while the edge cqi is parallel to the ^-direction, the coordinates of the rotated vertices are 
constrained by the following set of equations: 

xvi = Ri3 = sin sin < yvi = R23 = cos -0 sin = zvi = R33 = cos 6 > (2-41) 

XV2 = Rii = cos V' cos (j) — cos 6 sin cj) sin < yy^ = R21 — sin ip cos (p — cos 6 sin cj) cos 9 = (2.42) 

= = sin 6* sin > (2.43) 

XV3 = R12 cos ^p sin + cos 9 cos cj) sin ip = yy^ = R22 = — sin ip sin 9 + cos 9 cos </> cos ^/^ > (2.44) 

zy^ = R32 = — sin 9 cos (/> = (2-45) 

By solving this set of equations we find that the Euler angles ip, (p s-nd 9 that give rise to the configuration 
in Figure [5] are 

i) 9 = {n + and ip = [p +1)'^ for ?i=odd, p=even and < < ^ 

ii) 9 = (n + 1)1 and ^ = (p + 1)| for n=even, p=odd and ^ <9 <2-k 

It follows that the arrangement of edges under scrutiny has measure zero in S'0(3). 

By a similar method it can be shown that whenever an edge lies in a plaquette, or is parallel to a 
plaquette, one of the Euler angles will have to be equal to ^ for n odd or even. Therefore, that arrangement 
will have measure zero. This is not so for the general arrangement (number 1) delineated above. However, 
for any such arrangement, the values for the orientation factor and subsequently the geometric factor G^y 
will always be the same and, in zeroth-order, will not lead to any changes of the expectation value of the 
volume operator. 

Hence, the only cases of interest — i.e. the cases with measure different from zero — will not lead to a 
rotational dependence of the expectation value of the volume operator in zeroth-order. This should not 
come as a surprise since the geometry of a regular 6-valent graph is such that to each edge there corresponds 
a co-linear one. Thus whenever the term Cg for edge e changes from -1 to 1, the term cr^/ of the co-linear 
edge undergoes the inverse transformation. As a consequence there will always be the same number of 
detg g/ g'/ (cr) = —4 and detg g/ g// (a) = 4, although the triplets involved will be different in each case. It 
follows that the overall value of the geometric factor G-^y remains constant. 

A similar reasoning holds for the 8-valent graph since here too each edge has a corresponding co-linear 
edge. Therefore, there will always be an equal number of o"g = 1 and = —1. This implies that, as in the 
case for 6-valent graph, when no edge lies on a plaquette, the value of the expectation value of the volume 
operator for each 8-valent vertex will be rotationally invariant. On the other hand, the orientation of edges 
in an 8-valent graph in which one or more edges lie in a plaquette, or an edge is parallel to a given direction, 
have measure zero in 50(3), as was the case for the 6-valent graph. However, as previously stated, it is 
precisely such cases that lead to a change in the value of the geometric factor G.yy. 

For the 4-valent case the situation is somewhat different since there are no co-planar edges. Those 
arrangements of edges with respect to the stacks of plaquettes that cause drastic changes in the values of 
the orientation factor are the following: 

1. No edge lies in any plaquette. In this case we obtain | detg g/ g"(o")| = 4 for all linearly-independent 
triplets. 

2. Each edge lies in a given plaquette. This gives | detg g' g"(o")| = 2 for all linearly-independent triplets. 

3. One edge is aligned with a given plaquette, one edge lies in a given plaquette, and the remaining edges 
do not lie in — and are not aligned to — any plaquette. In this case we obtain | detg g/ g" (it)| = 1 for two 
triplets , I detg g/ ^t, ((t)| = 2 for one triplet, and | detg g/ g// (o")| = 4 for the remaining triplet. 

Similar calculations to those for the 6-valent graph then show that the cases 1 and 2 above have measure 
zero in 50(3). 
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In summary, the discussion above shows that for all 4-, 6- and 8-valent graphs, those orientations of the 
edges with respect to the stacks that cause a drastic change in the orientation factor have measure zero in 
50(3). Therefore, up to measure zero in SO (3), the geometric factor G^y for these graphs is rotationally 
invariant. 



2.4.1 Computation of the geometric factor for 4-, 6- and 8-vaIent graphs 

In this Section we will compute the geometric factor G^^y for the 4-, 6- and 8-valent graphs. We recall from 
equation (j2.33p that the expression for the geometric factor is 



G 



7,1) • 



1 

48 



e(e,e',e") det^^{a) 



e,e ,e 



(2.46) 



1 



^ e(ei,ej,efc)^det (a) 

l<i<j<k<N e,,ej,efc 



where A'" is the valence of the vertex and det(,.^(,j,e^i(^) = ^/Ji^'^ei^ej^el- what follows we will calculate 
G^^v for the 4-, 6- and 8-valent graphs respectively. In particular (for each valence) we will analyse each of 
the cases discussed in the previous Section which lead to different values of orientation factor. Any sub-case 
of these cases will lead to the same geometric factor. 



4-valent graph: We now compute the geometric factor for the 4-valent vertex for different embeddings 
of the graph in the stack of surfaces. 

1. The most general situation is one in which none of the edges is aligned to, or lies in, a given plaquette. 
Thus, for example, consider the situation in which the edges ei, 62 63 and 64 are in the octants A, C, 
H and F respectively (see Figure [7]). Such a combination has a non-zero measure in SO (3). 

The values for dete,,ej,efc (f) relative to this case are given by 



det (a) = det (a) = -4 

61.62.63 61,63,64 

det (fj) = det (cr) = 4 

61.62.64 62,63,64 



Inserting these values in (12.460 gives 

G'y^l] 



l<i<j<fc<4 



ei,e,-,efc) det [a] 



€-1 ,6^,6;^ 



1 



4e(ei, 62, 63) - 4e(ei, 63, 64) + 4e(ei, 62, 64) + 4e(e2, 63, 64)) 



-16 



(2.47) 



(2.48) 



2. If instead we consider the case in which each of the edges lies in a plaquette as, for example, is depicted 
in Figure [8l then the value for the geometric factor is 



G-y r) 



^ e{ei,ej,ek) det {a 



l<i<j<k<4. 



(2.49) 



2e(ei, 62, 63) - 2e(6i, 63, 64) + 2e(6i, 62, 64) + 2e(62, 63, 64)) 
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Figure 7: General 4-valent vertex 




3. For the situation in which one edge Ues in a plaquette and another edge is ahgned with a plaquette in 
another direction (Figure [9]), we obtain 



l<i<j<k<4 



ei,e,-,efc) det {a] 



(2.50) 



le(ei, 62, es) - 4e(ei, 63, 64) + le(ei, 62, 64) + 2e(e2, 63, 64)) 
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Figure 9: General 4-valent vertex 



However, we proved above that the embeddings of the vertex with respect to the stack depicted in cases 2) 
and 3) have measure zero in 50(3). 



6-valent graph: We now compute the geometric factor for the 6-valent vertex in the cases 1 to 4 described 
in the previous Section and which lead to different values of the signature factor. 

1. We start with the most general embedding of a 6-valent vertex with respect to the stacks. For 
example, consider the case in which the edges ei, 62 63, 64 65 and eg are in the octants A H E B C 
and G respectively. We then obtain the following value for the geometric factor: 



G. 



- E 



e(ei,e,-,efc) det [a] 



l<i<j<k<4: 



(2.51) 



;(4e(ei,e2,e3) + 4e(ei, 63, 64) + 4e(ei, 64, 65) - 4e(ei, 62, 65) - 4e(e2, 63, 64) + 4e(e2, 65, ee) 



-4e(e3, 64, ee) - 4e(e4, 65, eg)) 



-4x8 



2. For the geometric factor when only one edge and its co-planar edge lie in a plaquette (Figure [3]) we 
obtain: G^^^ = 2. 
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3. For the case in which two edges and their co-planar edge he in two different plaquettes in two different 
directions, while the remaining edge and its co-planar edge are aligned with the plaquette in the third 
direction (Figure [5]) , we obtain G^^^ = \/2. 

4. For the case in which all the edges are aligned with the stacks (Figure [6]) we obtain G^^y = 1, since in 



that case det 



1 for all linearly independent triplets ej, e,-, e^. 



However, we have proved above that cases 2), 3) and 4) have measure zero in S0{3). 



8-valent graph: We now compute the geometric factor for the 8-valent vertex for different embeddings 
of the graph with respect to the stack of surfaces. 

1. In the most general case, none of the edges lie in, or are aligned to, a given plaquette: for example, 
when the edges ei, 62 63, 64 65, cq 67 and eg are in the octants B, C, A, D, H, E, G and F respectively. 
This leads to the following result 



G. 



^ e(ej,ej,efc)^det {a) 

l<«<j<fc<4 " 



(2.52) 



-(4e(ei,e2,e3) + 4e(ei,e2,e4) - 4e(ei,e2,e8) - 4e(ei,e2,e7) - 4e(ei,e3,e4) + 4e(ei, 63,66) 
o 

-4e(ei, 63,63) + 4e(ei,64,66) - 4e(ei,e4,67) + 4e(6i,66,67) + 4e(6i,66,68) - 4e(6i,67,68) 



-4e(e2, 63, 64) - 4e(e2, 63, 65) 



-4 X 32 



(2.53) 



2. A more restricted case is when one edge and its co-planar edge are aligned with a plaquette in a given, 
different direction, while the remaining three edges and their co-planar edge lie in a given plaquette. 
Here we obtain G-^y = \/5- 

3. A special case is when each edge lies in a given plaquette: this gives G-yy = 2\/2. 

Similarly to the 4- and 6-valent vertex above, arrangement 2) and 3) have measure zero in SO{2>). 

From the discussion above of the geometric factor we can already deduce that, ignoring off-diagonal 
entries of the edge metric A, the expectation value of the volume operator gives the correct semiclassical 
value only for combinations of edges that have measure zero in 5*0(3). 

In fact, in zeroth-order in | the expectation value of the volume operator is given by 



J) 



dei{E){v)\ [dQi{dX{s)/ds)] 



X{s)=v 



1 
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^ e(6,6',6") det^/a) 

• • g g' g" 

er\e'r\e"=v 



(2.54) 



where {^)'^ ^J\dei{E){v)\\[dei{dX{s) / ds)\x[s)=v\ approximates the classical volume Vy{E) as determined 
by E of an embedded cube with parameter volume {al/b)^. It is straightforward to see that the correct 
semiclassical behaviour is attained for G^y = 1. 

The fact that the correct semiclassical behaviour of the volume operator is attained only for cases in 
which the graph is aligned to the plaquettation (6-valent case), or each edge lies in a given plaquette (the 
4-valent case), seems rather puzzling since, both cases, have measure zero in SO{2>). This makes one question 
the prima facie validity of utilizing the area coherent states to compute the expectation value of the volume 
operator. However, it is interesting to note that case 4) of the 6-valent graph is precisely what one gets 
when constructing such a graph as the dual of a cubical cell complex. 

We will now proceed to compute the higher, |-order dependence of the expectation value of the volume 
operator for 4-, 6- and 8-valent graphs respectively. 
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3 The Higher, |-order Dependence of the Expectation Value of the Vol- 
ume Operator 

In this Section we compute the expectation value of the volume operator for the 4-, 6-, and 8-valent graphs. 
The construction of these graphs was discussed in [23] in terms of regular simplicial, cubical and octahedronal 
cell complexes respectively. 

The following Section is subdivided into four parts. In the first we explain the general method to be 
applied in the subsequent Sections. In the second, third and fourth parts we apply this method to our 4-, 
6- and 8-valent graphs respectively. Each of these subsections is itself subdivided into three parts: in the 
first, the stack family and the cubulation that defines the platonic-body cell complex dual to the graph are 
aligned (see [23j); in the second we study the effect of a rotation; and in the third we study the effect of a 
translation. 

3.1 Initial preparations 

As a first step towards computing the expectation value of the volume operator, we must calculate the 
values of the quantities tl and ij^, defined in [23], which indicate the number of surfaces, s^j, that the edge 

e intersects, and the number of surfaces, s^j, which are intersected by both edges e and e'. Both these 
quantities depend explicitly on how the graph is embedded in the stack family (see |23j). In fact, the 
conditions for two or more edges to intersect a common surface are the following: 

1) Two edges Cj and ej intersect the same plaquette, s^^, iff < 0i, < ^ or ^ < (pi,(pj < n 

2) Two edges Cj and ej intersect the same plaquette, s^^, iS < 6i,9j < ^ oi ^ < di,9j < ^ 

3) Two edges and ej intersect the same plaquette, s^^, iff < ^i, 0j < vr or tt < 9i,9j < 2ir 



4) If we have equalities in any of the above conditions, such that the angles of each of the two edges 
correspond to a different limiting case, we obtain t^^^. = {ti G R\Sl H Cfc / 0, = = 0. 

We can also have situations in which two or more edges intersect a common plaquette in more than one 
stack. The conditions for such occurrences are the following: 

a) Given condition (1), two edges Cj and ej will intersect more than one z-stack iff \9i\ + \9j\ < 7r/2 and 
such that nf < 9i,9j < (n + l)f , where n = {1,2,3,4,5,6,7,8}. 

b) Given condition (2), two edges and ej will intersect more than one x-stack iff condition (1) above is 
satisfied and n| < 9i,9j < {n + 1)|, where n = {1,2,3,4,5,6,7,8}. 

c) Given condition (3), two edges and ej will intersect more than one y-stack iff condition (1) above is 
satisfied and n| < 9i,9j < (n + 1)|, where n = {1,2,3,4,5,6,7,8}. 

The conditions above imply that rotating the graph will change the values of the t^.^. and also the 
number of the ti „ that are non-zero. 

We will now briefly explain, with the aid of an easy example, the strategy we use to compute the terms 

tie- 

I ^ ^ ^ that are used in the calculations of the expectation value of the volume operator for the 4-, 6- and 

8-valent graphs. To this end, consider an edge, e € 7, of a generic graph whose length is given by 6. This 
edge will intersect the stacks of plaquettes in each direction a certain number of times. In particular, given 
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a length I of a plaquette, each edge will have n intersections with the stacks of any given direction, where 
n is identified with the Gauss bracket [y] and Cj is proportional to 5, where for i G {x, y, z} are the 
coordinates of the edge. 

For example, in the two-dimensional case of Figure [TOl the values of n in any given direction for vertex 
Vi (or equivalently the edge eo,i of length 5), whose coordinates are Vi = {6cos{a), 5sin{a)), would be 

vA 



1 21 X 

Figure 10: Intersections of an edge with the plaquettes in 2-dimensions 



Concomitantly, the expectation 



ii'x — [ I \ anu iiy — [ ^ J. 

The values n^, i G {x, y, 2}, depend on both the angle a and the ratio 
value of the volume operator will also depend on such parameters. 

The rotational dependence will be dealt with later. In the present Section we will focus on the 
dependence. We need to consider three different sub-cases: 



1. 
2. 
3. 



> 1 



< 1 



and determine which of the them leads to consistent solutions. 

However, to obtain an expansion of \[A we need to perform a Taylor series. The condition for applying 
such an expansion is that ||^ — 1|| < 1. From the expression for the (square) matrix A (see (j2.26p ) it is clear 
that the condition above is satisfied iff m(m — 1)-^= < 1 where m is the dimension of the matrix. As we will 



show. 



C X ^ where C=constant and /' < /; thus the condition m{m — 1) 



< 1 becomes I « 6e 



i.e. we need to choose the parquet to be much finer than the edge length (see Section [2]2]). If this requirement 
is satisfied, then we can perform a Taylor expansion of \fA obtaining \fA = 1-|-^(A—1)-|-|(A— 1)^-1-0(^—1)^. 
Actually, we are only interested in first-order terms, and so we shall only consider the approximation 
\J~A. ~ 1 -|- ^(j4 — 1) whose inverse, in first-order, is simply (\/]4)~^ — 1 — ^(^ — 1). Since the parquet length 
must be much finer than the edge length, in the remainder of the paper we will consider only case (1) and 
analyse whether it gives the correct semiclassical limit. 
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The first step in the calculation is to determine the range of allowed positions for each vertex, Vi, of the 
graph with respect to the plaquette. Since the graphs we consider are regular, determining the position of 
one vertex suffices to derive the positions of the remaining vertices in the graph. 

As an explanatory example let us consider a regular 4-valent graph 7 whose vertex Vq coincides with 
the point (0,0,0) of the plaquettation and whose vertex V2 (equivalently the edge 69,2) has coordinates 
-^). It follows that the range of allowed positions of V2 is from (nl, nl, nl) to (nl+1, nl+1, nl+1) 

where Ux = n. 



n 



[-^], as depicted in Figure [TTl 




(nl+1, nl+l,nl+l) 



(nl,nl,nl) 



Figure 11: Allowed positions of a vertex 



It is straight-forward to understand that different positions of V2 will determine different values of te^e for 
any two edges and Cj incident at V2. A detailed analysis shows that the terms te,ej differ according to 
which of the following conditions is satisfied: 

I) IxyJ > \nl + || 

II) IxyJ < |n/ + || 

III) Ixyjl = |n/ + || 

Similar conditions apply for all vertices in 7. 

Since the position of V2 will determine the positions of all other vertices, it is possible to establish which 
positions of V2 will lead to different values of the terms .'''"^ for all edges of all vertices of the graph 7. 

Such positions of V2 for a regular 4-valent graph are: 



a) 


\nl\ < \xv2 


< \nl + || 


b) 


\nl + || 


< 




< \nl + i\ 


c) 


\nl + {\ 


< 


XV2\ 


< \nl + || 


d) 


\nl + || 


< 




< \nl + || 


e) 


\nl + i\ 


< 




< + 


f) 




< 


\XV2 


< \nl + 21 
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For each such condition it is possible to derive the respective conditions for both the y- and the ^-coordinates 
in the three-dimensional case. It turns out that similar relations hold for the 6- and 8-valent graphs as well. 

To explicitly compute the terms t(.^ej we must choose one of the above conditions (a /), each of 
which will lead to different values for each teiCj- However, the computation procedures are the same. In the 
calculations of Sections 13.21 we will choose case (a). 

To describe the method for computing the values of teiCj ; we go back to a very simple example in 
two dimensions. We will then give the general outline of how this calculation can be generalised to the 
3-dimensional case. 

Let us consider Figure \T2\ where we chose x^^ > nl + ^. For simplicity we are assuming that the vertex 



y 



\ 
















V 






) tl 




1 






e . 














\ — 



nl nl+1 



X 

Figure 12: Example in two dimensions 



is symmetric with respect to the axis, i.e., the angles, (f), made by the two edges with respect to the y-axis 
are the same. 

We now want to compute the values of teiSj = teiSj + t^-Sj + ^e^ej where for each t^-g., k = {x, y, z} we 
have p. = i^. n i^.. 

As a first step we compute for each edge, e,, the value of i|. in the y-direction, obtaining 

= (nl + I — x^j) cot a and = (xt,^ — nl) cot a (3-1) 

Since the two edges commonly intersect only one y stack, in order to define the value of tl^gj = ili H we 
need to establish which of the two terms or tg. is the smallest. Thus, for example, 

Xy^ — nl > nl + 1 — Xy^ iff x > nl + - (3.2) 
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Since we have chosen x„j > n/ + 1 it follows that tl, < tf^ which implies that t|,ej = tij = {nl + l — x^^ ) cot a. 
As can be seen from Figure W2\ there are no intersections in the x stacks, therefore we obtain 



We now want to determine the values for 

sjEj (nl+l— x^-^) cot a 



where, in this situation, = Te +'^< 



(3.3) 

5sina+5cos a = 



te,] therefore, ,/ ] — , — . 

~\J '^i \/((5 sin a+(5 COS a)^ 

This calculation is very simple since the intersection of the two edges occurs only in one y stack. But it 
could well be the case that the angle between two edges is such that they intersect more than one stack in 
a given direction. For example, consider Figure [T3l always in two dimensions. In this case we would have 





X 



Figure 13: The second example in two dimensions 



Since in analysing the expectation value for the volume operator we will be considering graphs formed 
by regular 4-, 6- and 8-valent lattice, it turns out that the angles Oi — the angle formed by the projection on 
the edge on the x-y-plane and the x-axis — and the angle, (/>j, with respect to the z-axis for any edge, are 
such that two or more edges can only commonly intersect at most one plaquette in a given direction. 

When generalising the procedure described above for calculating the values of .J'"^ to the 3-dimensional 

case, some extra care is needed. In fact, consider Figure [HI It is clear that the values for t^. can be computed 
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with respect to both the x- and the y-coordinates as follows 

--axZl = c X cot 4)^^ = cx Zl (3.4) 

' ' smt' ' 

iff the edge points in the negative x direction 
iff the edge points in the positive x direction 

iff the edge points in the negative y direction 
iff the edge points in the positive y direction 

The term ^t^. in these equations represents the value of t^. as computed with respect to the x-coordinate, 
while ^tg. is the value of tg. as computed with respect to the y-coordinate. The non-uniqueness of the 
computation of the values t^. implies that there is an extra difficulty in the three-dimensional case. We will 
illustrate this with the aid of an example. 

Consider the edge in figure [T^l The value of t^^ can be computed with respect to both the x and 
the y coordinate, thus obtaining ^t^. or ^ig, respectively. However, it is clear from the diagram that the 
intersection of the edge with the stack of plaquettes in the z direction containing the vertex Vj is given 
by ^tg . . On the other hand ^tg . defines the intersection of the edge with the stacks of plaquettes in the z 
direction containing the vertex plus the stack in the z direction delimited, in the x direction, by the values 
nl + I and nl + al. 

This example shows that, given the values ^t^. and ^tg,, the intersection of the edge Ci with the stacks of 
plaquettes in the z direction which contain the vertex Vj is given by the smallest term, i.e., t^. =^ tg, tg. . 
It follows that, given two edges and ej, in order to find tg^ g^. := tg. Htg^. we first need to establish whether 
tg. =^ tg^ or tg. =y tg. and similarly for the edge ej. Once the value of the terms tg. and tg^. is determined, 
we can proceed as for the two-dimensional case and identify tg.g^. with the smallest t^, i.e. tg, g^. := tg. n tg^. . 

For intersections in the x and y stacks the procedure for computing the values of tg.g . is essentially the 
same. However, the formulae for the values of the individual terms, ■'tg,, are different. Specifically, for the 



''tl = a X cot (j)^— 
cos 



where 

a 



nxy.l + I - Xy. 



and 

I Uyy, l + l- yv, 
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(nl+al, ml+1, pi) 



(nl+1, ml, pi) 



Figure 15: Computation for the terms t^. in 3-dimensions 



x-direction we have: 

^i^^ = d X tan cos = d X F^^ ^tl^ = cx cotO = cx (3.5) 

where 

— i^zy, I iff the edge points upwards 
^zy. I + 1 — zy. iff the edge points downwards 
and c is defined as above. For the y-direction we have 

^t^^ = d X tan sin 6* = d X Fl = a x tan 6^ = a x Tl (3.6) 

where a and d are defined as above. 



d 



When computing the values of .^'''^ in three dimensions, as for the two-dimensional case we need to 

V "^i ^3 

compute the values for fg. , which in this case are simply te^ = t^. + te- + tg. = cos (90 — (pej cosOei + 
(5ei cos(90 — (/>£.) sin + cos (pi. The explicit values of the terms ."/^l obtained for the 4-, 6-, and 
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8-valent graph which satisfies condition (a) above, namely \nl\ < jxval ^ \nl + ^|) for the 4-valent graph and 
an equivalent condition for the 6- and 8-valent graphs are given in the Appendix. 

Since for all 4-, 6- and 8-valent graphs we are dealing with symmetric lattices, after a certain number 
of vertices the values for the terms ."''"^ will repeat, i.e., there will be a periodicity in the values of the 

terms ^==. Therefore, in computing these values we need only consider those vertices which comprise 

the periodicity cell, i.e., those vertices for which the values of the term ."''^^ cannot be obtained through 

V '^i "^3 

symmetry arguments. As we will see later, this periodicity is different for graphs of different valency. 

We now proceed to compute the expectation value of the volume operator for the 4-, 6- and 8-valent 
cases, utilising the values of the terms given in the Appendix. 



3.2 Analysis of the expectation value of the volume operator for a 4-valent graph 

In this Section we will compute the expectation value of the volume operator as applied to a 4-valent graph. 
We will first take into consideration the non-rotated graph. In establishing rotational and translational 
dependence of the expectation value we will perform both a rotation by arbitrary Euler angles and a 
translation and, then, recalculate the expectation value. We will see that the contributions that come from 
the terms — t==, which comprise the off-diagonal elements of the matrix ^/A , are not trivial, thereby 

producing a strong rotational and translational dependence in the expectation value of the volume operator 
in higher order in |. 
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3.2.1 Expectation value of the volume operator for a 4-valent graph 

To calculate the expectation value of the volume operator we will consider a 4-valent graph constructed 
from the simplicial cell complex as discussed in [23]. We choose the vertex Vq to be Vq = (0,0,0), and the 
angles (j)e = ^°^~^("^) ~ ^'-'^ all e G 7 such that we obtain the configuration depicted in picture 

[T6l The periodicity cell for a 4-valent graph contains four vertices, including Vq. The coordinates of the 
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1 

Figure 16: General vertex of a 4-valent graph 



remaining three vertices are V, = (^,^,^); Vs = (^,3^,^); ^13 = (0,2^,2^). 

It follows that the edges eo,i, eo,2i ^0,3 and eo,4 lie in the octants G, B, E and D respectively. This 
implies that the geometric factor for the vertex Vq will be G^y = y/2. Because of the geometry of a regular 
4-valent graph, it turns out that all the vertices comprising the periodicity cell, V2, and V13 will have 
G-iyi = \/2. 

The following table gives the values obtained for the terms "''"^ for the 4-valent graph that satisfies 

condition (a) as defined in the previous Section, namely |n/| < [xval < |n/ + ||. 

It should be noted that, because of the geometry of the 4-valent graph, the terms Tei'^^ ^Fe^'^^ ^Ze^'^^ in 
equations 13.41 13.51 and 13.61 are all equal to 1 for each edge ej : 



^0 


six terms 




= 








V2 


six terms 












six terms 












two terms 




( 2(5 









Here, 5 is the length of the edge e. In order to apply equation (|2.20p . we first need to determine the values 
of the term det^ g/^// (\/]4 ) for each triplet of linearly-independent edges e, e',e". Using the fact that, in 
first-order approximation, (\/A)~^ = 1 — \{A — 1), the explicit expression for (ys/~A)~^ for the 4-valent graph 
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under consideration is 
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where a 



V3 



-nl) 



—-7= and the terms e > 



and are computed using the techniques defined 



in the previous Section. 

Now that we have an expression for the inverse of the matrix ^/A we can compute the expectation value 
of the volume operator for each of the four vertices in the periodicity cell and, then, sum their contributions. 

We start with the vertex Vq. First consider the sub-matrix of the matrix {■\/~A)~^ formed by all the edges 
incident at Vq. This is 

60,2 



60,1 
60,2 
60,11 
60,12 



60,1 



60,11 60^12 \ 





1 










1 










1 



Because of the geometry of the 4-valent graph, at each vertex there are four triplets of linearly-independent 
edges. Keeping this in mind and computing the determinant of the matrices formed by each such set of 
triplets, we obtain the following expression for the expectation value of the volume operator at Vq: 



5^/^./det iEHu)) IGldetf 



By a similar procedure for vertices V2 and V3 we obtain 



det [E^Au)) 



16 1 



3a^ a-^ 



4 4 

In both cases, the sub-matrix of \/~A ^ we consider is 

/ 1 



(3.7) 



(3.8) 



1 , 



1, 



— |q; 1 —^ct 

1 . 



2 

-ha 



? 
2 

1 y 



a 



and then we compute the determinant of all the sub-matrices formed by linearly-independent triplets of 
edges. 
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For the vertex V4 we obtain 



where we have used the sub-matrix 

/ 1 



Summing up these contributions gives 



16(1 - 



det 



-a 



d{s, u^,u'^] 
\ 



(3.9) 
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Vr = 6' 



det 



+ J\det {Enu))\(\\l 



3a2 



1 J| det {E^{u)) I + J| det {E^{u)) | ( | 1 



a3 A I 



4 4 

To first-order approximation we obtain 



+ J|det {E-iu))\(\l 



a 



Vr = S'^lVl J \ det {EVu)) 



det 



d{s, u^,u'^^ 



1 /3a^ 



2 V 16 



0{3) 



(3.10) 



(3.11) 



It should be noted that, although the term det(£'j'(u)) is vertex dependent, we can safely assume that, 
to first-order in 1/6, the values will be the same for each vertex within each periodicity cell that involves 
only an order of four vertices. Thus this term can be factored out from the equation. This first-order 
approximation will be used throughout the rest of this paper. As mentioned previously this is justified since 
we choose | << 1. It follows that the terms oc which are much smaller than one (see Section [2. 2p . 

The term proportional to in the equation above represents the I / (5-correction for the expectation value 
of the volume operator for a given region R. As in [23], for a general 4-valent graph, even in the zeroth-order 
approximation the expectation value for the volume of a given region R does not coincide with the classical 
value for the volume of that region. Notably, there is no linear correction in l/6\ 

3.2.2 Expectation value of the volume operator for a rotated 4-valent graph 

We now analyse how the results of the calculations above depend on how the graph is embedded in R^. 
Here we will consider rotational invariance; translational invariance is discussed in the following subsection. 

To analyse the rotational dependence of the expectation value of the volume operator, we will perform 
an Euler rotation of the graph with respect to some arbitrary Euler angles /3, ijj, a and then repeat the 
calculation. The transformation matrix is 

cos (j) cos ip — sin (/) cos 9 sin ip cos ip sin <j) + cos 9 cos <j) sin ip sin ip sin 9 
i? = I — sin tp cos (j) — cos 9 sin (p cos ip — sin tp sin (p + cos 9 cos (p sin tp cos ip sin 9 
sin 9 sin (p — sin 9 cos (p cos 9 

The coordinates of the rotated vertices are then given hy V- = R - Vi 

RiiXv, + RiiVvi + RizZv^ 
Rsi^v, + R?,2yvi + -Rss^vj 
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Applying this transformation matrix to the 4-valent graph we obtain the following new coordinates for 
the vertices: 



s s s 

Vi = ((-i2ii+i?i2 + i?13)^,(-i?21+i?22 + i?23)^,(-i?31+i?32 + i?33)^) (3.12) 
^8 = ((i?ll+3i?i2 + i?13)^,(i?21+3i?22+i?23)^,(i?31+3i?32+i?33)^) (3.13) 
1^/3 = (2(i?i2-i?l3)^,2(ii22-ii23)^,2(i232-i233)^) (3.14) 

The new angles between the rotated edges and the x,y,2:-axes can now easily be computed using elementary 
trigonometry. 

As an explanatory example let us consider the edge eo,2. To find the angles this edge has with respect 
to the axes we need first to compute the coordinates of the vector eo^2 starting at vertex Vq and ending at 
vertex ¥2- In this case, the coordinates of 69^2 coincide with the coordinates of the vertex ¥2'- 

s s s 

62,0 = [{Ru + R12 + ^13)^, (^21 + R22 + ^23)^, (^31 + ^32 + ^33)^) (3.15) 

If instead we considered the edge 625 we would get 

S 3 S 

625 = eo2 - eo5 = ( - {Ru + R12 + ^13)^, -(^21 + R22 + ^23)^, -(^31 + ^32 + ^33)^) (3.16) 

Once we have the coordinates for 62,0, the angle, (/>e2 q, it forms with respect to the z-coordinate is 



, , -ifV^T^\ , 1 f Vi-Ru + ^12 + Ris? + (-^21 + R22 + R2W , ,.,7, 

'--='^" (^^)=*^^ I (-i^3,+i^32 + i?33) ^ ^'-''^ 



The angle, 0^^^, between the projection of the vector on the x-y plane and the x-axis is given by 

0,, „ = tan-i = tan-i ( Z^21±^3l±^) (3.I8) 

In the same way we can obtain the angles for all the edges in our graph in terms of the elements of the 
transformation matrix. Thus the orientation of each of the edges of the graph will depend on the matrix 
elements of the transformation matrix, i.e., on the Euler angles that parameterise the rotation. 

In order to determine the rotational dependence of the expectation value of the volume operator, we 
have performed a case study in which the expectation values were computed for all possible orientations 
of the graphs that have non-zero measure in SO (3). Such possible orientations were described in Section 
12. 4[ To aid calculational simplicity, these sub-cases are defined in terms of possible ranges of values for the 
angles and (p for each edge in the graph, rather than on possible values for the Euler angles. 

In order to keep our results as general as possible, we performed our subdivisions so as to cover all 
possible situations. This is less tedious than it might seem since we are dealing with regular lattices and, 
therefore, once the angles for the edges of one vertex are fixed, we immediately know the orientation of the 
edges of all other vertices. 

Let us choose Vq as our reference vertex with respect to which the possible orientations of the edges are 
defined. The edges incident at Vq are eo,i , eo,2 , 69,3 , 60,4. In defining the orientation we use the convention 
that both < < 27r and < 9 < lir increases anti-clockwise. 

Once the rotational matrix has been applied, whatever the values of the Euler angles might be, we will 
end up in a situation in which two edges e,, ej point upwards, i.e., — ^ < (pei, (pej < f j and the remaining two 
edges point downwards, i.e. ^ < (pe^., 4>ei < This is a consequence of the geometry of the 4-valent graph. 
We will call two edges pointing in the same up, or down, direction an 'up', or 'down', couple respectively. 
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Since we are considering only those edge orientation with measure non-zero in SO (3), the angles of the edges 
Cj, Cj of each up/down couple will satisfy the following conditions: |^ej = |90 + ^ej | and {(pal = \4'ej — 54.75| 
Given a particular choice of up and down couple we have to specify in which octant (see Figure H]) 
each edge lies. This is required since different octants induce different values for the geometric factor G^y. 
The angles 6(,^ and (j)^,^ required for an edge ej to lie in each of the octants are listed in Table [1] where, 
again, we use the convention that Q < (pei < 27r and < < 27r, with both angles increasing in an 
anticlockwise direction. However, because of the geometry of a 4-valent graph, the allowed angle-ranges 





A 


B 


c 


D 




^ < 0,. < 27r 


< < f 


< (t>e, < f 






< fe, < 1 


f < 6)^, < TT 




f - < ^e. < 2lV 




E 


F 


G 


H 




TT < 0,, < if 


f < 0^, < TT 


f < (?!)e, < TT 


TT < (jie, < f 




< fe, < f 


f < fe, < TT 


TT < fe, < f 





Table 1: Angle-ranges for each octant 



have to be restricted to those listed in Table [21 
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D 




^ <(l>e, < 27r-sm-^(i 


sm-^(i < 0.. < f 




f < (jie, < 27r-sm-^(i 




< ^e. < f 











E 


F 


G 


H 




TT + sin-'d) <(t>e, < ^ 


f < (/>e, < TT -Sin--'(i) 


f < <;/>ei < TT - sin--'(i) 


7r + sin--'(| < (jie, < ^ 




< e., < 1 


f < 6»e, < TT 




^ < 6*,. < 27r 



Table 2: 4-valent graph angle-ranges for each octant 



Our calculations show that for all possible sub-cases of angle arrangements defined in Table [2l the 
expectation value for the volume operator is rotational invariant only at the zeroth-ordei|£|, while higher- 
order terms are rotationally dependent. Therefore, in what follows, we will not compute the expectation 
value for the volume operator as computed for each possible orientation of the graph. Instead we will choose 
a particular sub-case of Table [2] and compute the expectation value for such a sub-case. Specifically we will 
choose the case in which the arrangement of edges, incident at the vertex Vq, after a rotation, is given by 
the following ranges: 

vr Svr _i 1 

< 6*60,2 < 2 Y ^ < ~ 3 

vr Svr 

- < Oeo,! < vr Y < <Aeo,i < vr - sm (-) 



Svr „ 1,1. Svr , , 

— < eeo,3 < 27r vr + sin \-) < 0^0,3 < — (3-19) 



This implies that the edges eo,2) ^0,3, eo,4 and eo,i lie in the octants A, H, C and F. Prom the geometry of 
the 4-valent lattice, the angles of the edges incident at all the other vertices follow. 

There is a vast range of Euler angles for which the case above is obtained but, for the sake of brevity, 
we will not list them here. What is important, though, is that such case has a non-zero measure in 50(3). 



^This is a consequence of the fact that the geometric factors G-y^Vi for each of the sub-cases in Table [2] will be the same (see 
Section I^XTj) 
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It should be noted that different combinations of angles within the angle ranges in (|3.19|) lead to dif- 
ferent outcomes for the expectation value of the volume operator, since they lead to different values of the 
terms — =|=. However, in zeroth-order, the expectation value of the volume operator will be the same 

irrespectively of which angles satisfying (I3.19P we decide to utilise. In fact, the rotational dependence of the 
expectation value of the volume operator, in the zeroth-order in |, is determined solely by the geometric 
factors G^y.. For the case which we are analysing (|3.19p . the values of G^yy. will be the same irrespectively 
of which sub-case of ()3.19p we analyse. On the other hand, the dependence of the expectation value of the 
volume operator on higher orders of 4 is determined by the terms .J'"^ and, therefore, will depend on the 
sub-cases we analyse. 

This discussion shows that for higher orders in | the expectation value of the volume operator is rota- 
tional dependent since, as stated above, for differing angle-ranges that lead to the same geometric factors, 
the values of the terms i",^"^, will differ. 



We will now compute the expectation value of the volume operator for the periodicity cell in the 4-valent 
graph for the case in which the angles of the edges incident at vertex Vq satisfy condition ()3.19p . The first 
step in order to compute the expectation value of the volume operator is to define the matrix , whose 

off-diagonal entries are the terms — ^===. This matrix is given in Table lA. 1.1 1 in the Appendix. Although 

different combinations of angles satisfying condition (j3.19p will lead to different values of the terms — ^==, 

however, any such combination will lead to the same non-zero entries of the matrix \/]4 . This means that 
the pairs of edges commonly intersecting a plaquette in a given direction will coincide for any combination 
of angles satisfying conditions (I3.19p . even though the number tei.cj of plaquettes they commonly intersect 

will differ in each case. Therefore in computing the matrix y/A ^ we will not determine the precise value of 
the individual entries, but we will leave them as general as possible. Their precise values can be computed 
once a specific combination of angles satisfying ()3.19p is chosen. 

Given the matrix ^/A we are then able to apply formula (I2.33P for computing the expectation value of 
the volume operator. As in the aligned case, we first compute the expectation value of the volume operator 
for each of the four vertices, and then sum their contributions. In what follows, the term .'""'7 is denoted 

by aei,ej = J^'^i ■ An explicit form for these terms can be found in Section 11.1 of the Appendix. 

V '^j 

The expectation value for the volume operator for the entire periodicity cell, up to first-order in | is: 



fl / BJC^ \ / I 



\ A, \-\ I ('— ry2 — (T^ — (T^ —0-2 — — it^ \\ /\rlpf(F'^ 

-T^y\±-r^\ "ei3, 18, 613,15 "ei3, 19, 613,15 "613,19,613,18 613,3,613,15 "613,3,613,18 613,3 ,613,19 7 1 Y I V^j 

+ 4^|1 + ^(-<,o,6,,5 - <,o,e.,6 - <,o,e,,r " <5,6.,6 " <5,e.,r " <6,e,,,)l ^l^^^l^i (^)) I ) + ^(3) 

(3.20) 

By performing a Taylor expansion for each of the roots present in the above formula, we can factor out the 
term ^\ det since, in the first-order approximation that we are considering, they turn out to be 

the same for each vertex. Such an approximation is justified by the analysis performed in Section (|2.2p . We 
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then obtain 



j\"'JJ\y-^ ' ^\ ^eg, 10, eg, 11 eg, 10, eg, 12 eg,io,eg,5 eg, n, eg, 12 eg, 11, eg, 5 eg, 12, eg, 5 
ei3,ig,ei3,i5 "ei3,i9,ei3,i5 '^ei3,i9 ,ei3,ig "ei3,3,ei3,i5 '^613,3 ,ei3,ig '^ei3,3,ei3,i9 {o.Zl) 

2_2_2_2_2_2\\ /o r)r)\ 

'^e2,o,e2,5 '^e2,o,e2,6 '^e2,o,e2,7 '^62,5,62,6 '^e2,5,e2,7 '^e2,6,e2,7'' J 

f(n-l) 

= 6-^^\det{E^iu))\2V2(l + \{- ^ 4) | det (^-^) | ) (3.23) 

The term ^i—YlIj=i-j^^^'ji) represents Z/(5-corrections. Each term aij is proportional to Cj for I' < I; C 
is a constant that depends on the Euler angles we chose. On the other hand the geometric factors G^y. 
for cases (I3.19P coincide with the geometric factors as computed for any of the sub-cases in Table [21 i.e., 
G^y^ = This implies that although for such cases the expectation value of the volume operator is 

rotational invariant in zeroth-order, nonetheless, it does not reproduce the correct semiclassical limit. 

For those embeddings whose measure is zero in 50(3), the geometric factor turns out to be different 
and, in zeroth-order in Ijb^ leads to a different value of the expectation value of the volume operator as 
computed for 4-valent graphs. 

3.2.3 Expectation value of the volume operator for a translated 4-valent graph 

In this Section we will analyse whether the expectation value of the volume operator for the 4-valent graph 
is translational invariant with respect to the plaquette. 

To perform this analysis we consider our original aligned graph and translate it by an arbitrary vector 
■p = {ex,ey,ez)- The new coordinates for the vertices are: 

Vq = {ex,ey,ez) with > ey > ez (3.24) 
Vi' = (- -I + + (3.25) 

V; = (-| + .„., + 3A,,, + J^) (3.26) 

Similarly to the analysis for rotational invariance, the computation of the expectation value of the volume 
operator can be divided into different sub-cases, each of which would lead to different outcomes. 

The first division is given by the choice of the signs and the relations between €x, ey and e^, i.e., whether 

they are positive or negative and whether one coordinate is bigger or equal to another. Each of these cases 

can be ultimately subdivided into sub-cases depending on the relation between the ratio b = A= and the 

v3 

coordinates of the translational vector. 

To carry out our calculations we choose the following: 

1) b> ex> ey> ez>Q 

2) iVfci - 1 1^/1 > 4/ - n{/ for ah \Vl\ > \Vi\ 

Altogether, such conditions will allow to determine both the sign of the coordinates for each of the vertices 
of the translated graph and also the magnitude relation between the coordinates of each vertex. 

However, it will transpire that our result is independent of which case we decide to use to perform the 
calculations. In fact, as we will see, in zeroth-order, the expectation value of the volume operator for a 
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4-valent graph is translation invariant up to combinations of measure zero in 50(3). However, for higher 
orders of approximation this wiU no longer be true. 

As a first step in our calculations we need to specify the allowed positions for each of the translated 
vertices. Due to the highly symmetrical structure of the 4-vaIcnt graph, in order to determine the allowed 
positions of each vertex, it suffices to find the allowed positions of one reference vertex. We choose such a 
reference vertex to be Vq, whose new coordinates are Vq = {^xi^y^^z)- 

The number of stacks intersected by the vector that represents vertex Vq in the x, y and z-directions 
arc, respectively, n = [y-], m = [y-] and p = [^] (where [] indicates the Gauss bracket). It follows that 
the allowed positions for vertex Vq are given by the following ranges of each coordinate: nl < €x < nl + I, 
ml < €y < ml + I and pi < ex < pi + I- 

It turns out that to carry out the calculations for the expectation value of the volume operator we have 
to restrict the value-range of the coordinates ex, Cy and e^. We choose nl < ex < nl + ^, ml < ez < ml + ^ 
and pi < ez < pi + ^ 

We will now compute the expectation value of the volume operator of the periodicity lattice of the 
4-valent graph. We will not give the details of all the calculations involved since they are quite lengthy. 
However, the method utilised is the same as for the non-translated case, namely, for each of the four vertices 
comprising the periodicity cell we consider the sub-matrix of ^/A labelled by the four edges intersecting 
at the vertex. For each of these sub-matrices, call them M, we compute the determinant of the four 3x3 
sub-matrices of M defined by the triplets of linearly-independent edges. We then sum up the contributions 
coming from each of the vertices. Similarly as for the aligned 4-valent graph we have 



The expression for the matrix -\/A ^ is 
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(3.28) 
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where Ai = {xy, - riil)^, Bi = {yy, - ml)^, Ci = {zy - Vil)^- 

Applying the method described above we compute the expectation value for the volume operator for one 
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periodicity cell to be 



Vr 



'1 



2x 



(3.29) 



det I ^J\A-B^-2C^-BoC^\ + ^ | det (^Ef{u)'^ \ ^ \A - 2Al - AIC2 - Cl\ + 

det {e-{u)) I ^14-2^2^-^236-13-0231 + ^1 det [e-{u))\^\A - ^ - 2C72 - A,Cl\^ (3.30) 



In first-order approximation we then obtain 

T., r N,/. 1 



=2W -6\ I det ^^(tx) I 4 - - LBo^ + 2Co' + BoC^ + 2Ai + A^Cs + + 2A{^ + AfgCia + C7i^3+ 



^1 + 2C7| + A8C| +0(4) 



det ' ^ 



(^(s, li^, li2) 



=2^25^1 det (^^.n^x))!!^! - ^ (i?o' + 2Co2 + 2^ + C| + 2A?3 + Cl^ + ^ 2C|) + o(3)^ 



det ' 



5(s, 
(3.31) 



where in the last equation we have only considered first-order contributions obtained by the usual Taylor 



series of the square root (see Section [T2]) . Thus, we were able to factor out the term y det(ii^^(u)). As it is 
evident, the corrections of second order in Ijb are not translationally invariant. 

3.3 Analysis of the expectation value of the volume operator for 6-valent graphs 

In this Section we will calculate the expectation value of the volume operator for a 6-valent graph. First 
we consider the non-rotated graph, then we will analyse the rotational and translational dependence of the 
expectation value by performing a rotation of the graph, followed by a translation of the graph. We will 
then recalculate the expectation value. 



3.3.1 Expectation value of the volume operator for a general 6-valent graph 

Similarly as for the 4-valent graph, we will analyse the case for which j > 0; the motivation for such a 
choice was given in Section [2. 2i For computational simplicity we will position the graph so that the (0, 0, 0) 
coordinates of the graph coincide with the (0, 0, 0) coordinates of the plaquette. We also need to align the 
graph in such a way that each vertex is symmetrical with respect to the axis. Therefore we will choose, for 
each vertex Vfc, the value (/>ej. ^ = 45 for all edges ek,i incident at Vk and 9e^ , = 45 for four edges, while the 
remaining two will have 9e^ ^ = 0. This edge orientation corresponds to the limiting case (2) described in 
Section [2X11 

As for the diamond lattice, we choose the vertex V2 as our reference vertex with respect to which we 
determine the allowed positions of all the remaining vertices of the graph. We also choose the allowed values 
of the j;-coordinate of V2 to be nl < < nl + ^, where, in this case, n = [^^]. Using the same method 
used in Section 3.1 we can compute all the terms '''"^ for the periodicity cell of the 6-valent graph that 
contains nine vertices. Given the geometry of the 6-valent graph we have the following values for the term 
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in equations 

\l iff^e, = [1 ifF^e, =0 



oo iff Be- = 



F^^ = 1 Vei G 7 (3.32) 

The coordinates of the vertices of the periodicity ceh are 



Vo = 


(0,0,0) 


V2 = 




Vs = 


((4-(4.-4.2.|: 


V4 = 




V5 = 


(0,-24,0) 


Vi3 = 


^ 2 ' 2' 2 ^ 


Vis = 


(26^ + 6^-16^-, 6^) 
^ 2 2 2 2 ^ 


Vi5 = 


(if -5,-24, Sif) 


= 


(4+4.-34,0) 



The values for the terms . ' I are given in the section IA.2.11 of the appendix 

The value for the expectation value of the volume operator for one periodicity cell consisting of nine 
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vertices is computed to be 
1 



Vr 



8 



det 



2 V8W|det fs^(n))| + 



23(a')2 a(a')2 3a(a') 



4^2 2Vldetii^>^ 



a'/3/3' a'^/3' 7(/3')2 A3(/3')2 



-^/2^U3/?' 



2^/2 



2^/2 



det (£;^(u))| + \/|8- 12/32 _ 8^/2/32 1 ^/ 1 det {E'^{u] 



+ V|8 - 2a2 - 9{a'Y - 2^/2a(a')2 - 12^2 _ §^^2 _ i2a'/3/3' - 9\/2a'/3/3' - 12(/3')2 - 4a(/3')2| J| det f£;«(n) 



+ 



I^^-Ii^-lM«Z_2«/3-2V2a/3-i^^-^-6/32-^ 
4 2 8^2 2 2^2 2^ 



2a'/3/3' -4(/3')2 



det (i^«(n)) I + ^/|8 - - a2/2 - - ^ - 3/32 - 2^2/32 - ^3/3/3' 



2,X^^(50 _ a 1313' 
2V2 2~ 



.^«M_3(/3')2_^OT 



4^2 



det (i?^(n)) I + J|8 - 3(4,)2 - (^^ig^li^ - 2^?, - _ 9(a')2 



9A,,(a')2 3(4,)2/3 (4,)2^ ,3V2^,./3-^ -10/32 -4^2/32- ^14^ -3a'/3/3'-3V2a'/3/3' 



4^/2 



^/2 



V2 



-2(/3')2|a/| det ( i?,»(n) ) | + J|8 - 2(A;,)2 _ _ _ ^A^^ _ 3V2^,,/3 



-I8/32 - 9^/2/32 - Mli^ _ _ 9(/3')2|y|det (s^(n) 



(3.33) 



Expanding the square root (see the analysis in Section I2.2p and considering first-order terms we obtain 



y^ = 25y|det {E-iu))\ 
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5(s, U^, ^2) , 

1^ _ 29(a')2 - 3V2A14/3 - 3^/2^l5/? - 2a/3 - 2V2a(3 - - 32^2/3^ - 

(3.34) 



Contrary to [23] we find that, to zeroth-order in 1/6, the expectation value of the volume operator for a 
6-valent graph does not have the correct semiclassical limit. On the other hand, if the graph is aligned to 
the orientation of the plaquette we do obtain the correct semiclassical value. However, this embedding has 
measure zero in 5*0(3). 



3.3.2 Expectation value of the volume operator for a rotated 6-valent graph 

We will now analyse the expectation value of the volume operator for a rotated 6-valent graph. As for the 
4-valent graph, different choices of Euler angles in the rotation give different values of teisj ■ Therefore, we 
will once again have to define sub-cases which are defined according to the possible ranges of values for the 
angles and Oe^ for each edge e^. Because of the geometry of the 6-valent lattice, we know that for any 
edge, e, of a given vertex there exists a co-linear edge, e', which intersects the same vertex. This implies that, 
given two co-linear edges e and e', we can define the angles of e (respectively e') in terms of e' (respectively 
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e) as follows: </>e = 180° — 0e' and 9^ = 180° — 9^'. Such relations reduce the number of cases that need to 
be analysed. 

We choose the vertex Vq as our reference vertex. The relations for the angles of the edges incident at Vq 

are: 





= 180° - 






= 180° - 


*^eo,i7 




= 180° - 




^eo,8 


= 180° - 


^eo,6 


*eo,i7 


= 180° - 


^eo,i 


^eo,2 


= 180° - 


^eo,7 



(3.35) 

These relations imply that the allowed values of the angles of the edges at a given vertex fall into one of the 
following groups: 

1. A given triplet of edges points upwards, i.e., their angle (j) lies between — | and |, and the triplet 
formed by their co-linear edges points downwards, i.e., their angle lies between ^ and This 
situation arises when none of the edges is aligned with one of the x, y, z-coordinates. However we 
have two distinct sub-cases which satisfy this arrangement of edges 

a. No edge lies in any plaquette. 

b. Only one edge and its co-planar lie in a given plaquette (Figure [3|). 

2. A given couple of edges points upwards i.e., their angle (p lies between — ^ and ^, and their co-linear 
edges point downwards i.e., their angle lies between ^ and This situation arises when one edge 
(and subsequently its collinear edge) is aligned with one of the coordinates axis and, subsequently, the 
remaining two edges and their co- linear lie in two different plaquettes in the same direction (Figure Ej) 

3. Only one edge points upwards, i.e., its angle cj) lies between — ^ and ^, and the co- linear edge points 
downwards i.e., its angle (p lies between | and This situation arises when all of the edges are 
aligned with the coordinate axis (Figure [6|) . 

A discussion of each of these cases and the respective value for the geometric factor was carried out in 
Section [2.4. 11 There it was shown that only case la above has non-zero measure in S0(3), therefore we will 
restrict our analysis to such a case. 

It is straight forward to see that case la can be divided into further sub-cases according to the values of 
the ^-angles and the relations between the i;A-angles of each of the up/down couples. In what follows, we will 
not give the results for all possible choices. Instead, we will choose a particular sub-case and perform the 
calculations for the expectation value of the volume operator with respect to this sub-case. As we will see, 
these calculations show that, up to embeddings of measure zero in S0{3), the semiclassical behavior of the 
volume operator, in zeroth-order, does not depend on how the graph is rotated: a fortiori, it is independent 
of the particular case we analysed. 

In order to carry out a proper comparison between the semiclassical behavior of the volume operator 
as applied to graphs of different valence, we will apply the same Euler transformations (i.e., with the same 
Euler angles) to each of the graphs we consider. Since we have not specified the values of the Euler angles, 
the only way to do this is to assume that after a rotation, those edges which had the same angles used in 
both the aligned 4-valent and 6-valent case will end up in the same octant. For example, consider Figures 
[T71 and [18] which depict both 6-valent and 4-valent vertices respectively. From such pictures it is easy to 
see that the edges eo,2 and eo,3 of the 4-valent graph have the same 9- and (/>-angles as the edges, eo,8) and, 
eo,i, of the 6-valent graph. Therefore, in the rotated case we will assume that eo,8) and, eo,i lie in the same 
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octants as 69,2 a-nd 69,3 respectively. It follows that the angle-ranges for the edges incident at vertex Vq for 
a 6-valent graph are: 



Svr 



2 < ^-0.2 < 4 



Svr 



< 



< 2tt 



IT 



TT 



'eo,i7 



< TT 



< 0eo,8 < 



< (/>eo.i7 < ^ 



vr 
2 

TT 



< 



^eo,7 < 4 



< 



'eo,6 



< 



Svr 

T 
Svr 



TT 

2< 



^eo,7 



< vr 



2 



vr < 



Svr 

T 

, Svr 



(S.S6) 



Such conditions of the angles implies that the edges eo,i, eo,2, eofi, eoj, eoj and eo,i7 lie in the octants H, 
D, G, F, A and i? respectively. 

Since we are considering a regular 6-valent lattice, the above ranges of angles induce a relation on all the 
other angle-ranges of the edges for each vertex in the graph. The coordinates of the rotated vertices are: 



^ — ■ \ V 
1 \ / 




1/ v\ 
/ o\ 




V 

1 , 


/ -> \- 











V 
6 



Figure 17: 6-valent vertex 




Figure 18: 4-valent vertex 
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vis 



vis 



VI, 



vL 



0,0,0 



{Rii + Ris)—^, {R21 + -^23)—^, (-R31 + 2i?33)— ^ 



(^ii( 



^)-^+i^23v/2)5,(i?3i(^ 



1, 



2 ' 2 

(^2^11 - i?i2 + V2Rn)6, (V2i?2i - i222 + V2R23)d, ( V2i?3i - i?32 + V2R33)6 



^ + R,,V2)6 



R12S, R225, R32S 
--^Rii - R12 + -^-^^3)^, {—-^R2i 



R22 + ^R23)S,i-^R 
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R32 + -^Rss)^ 



{{V2 + ^)Rn - ^Rl2 + ^Rl3)S, {{V2 + ^)i?21 - ^ii22 + ^R23)S, {{V2 + ^)R3l - ^R32 + ^i?33)'^ 



V2. 



V2 
2 

.\/2 + l 



1)R 



11 



i?12 + ^Rl3)S, ((^ - l)i?21 



R22 + ^R23)S, - 1)^31 - i?32 + ^i?33)5 



R 



11 



^^12)<^, (^^^21 - Ir22)5, (^^i?31 - ^ii32)5 



(3.37) 



Similarly, as was done for the 4-valeiit case, different choices of combination of angles satisfying conditions 



(j3.36p above will lead to different values for the terms 



. However, the couples of edges commonly 



-1 



intersecting a given stack will coincide for any such combination. This implies that the matrix will 
have the same entries for any sub-case of (|3.36p but their specific values will be different. 

Moreover, the geometric factor G-yy of any sub-case of (j3.36p will coincide. It follows that any com- 
bination of angles that satisfies conditions ()3.36p will lead to the same value in zeroth-order in | of the 
expectation value for the volume operator. Therefore, as was done for the 4-valent case, in order to compute 
the expectation value for the volume operator, we will not specify a particular sub-case of (I3.36p . but leave 
the result as general as possible. 
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Given conditions 13.361 the expectation value of the volume for the periodicity cell is 



1 



"^■\/l^~^9i ^^63,83 /'63,23 /'lla.es ^113,83 ^113,133 ^113,23 ^173,63 ^173,83 /'l73,133 Pn3,23 



-/'i33,83 -/'i33,23 jl ^ \l\det(^E^(u) 



Pl4i3,21i3 /'9i3,14i3 /'3i3,21i3 /'3i3,9i3 ^'4i3,9i3 ^'4i3,21i3 ^'4i3,3i3 ^'4i3,14i3 /'5i3,9i3 



-P5i3,21i3 - P5i3,3i3 - P5i3,14i3 I >< A/lc^etLElH 



+ 7: 



2V ^02,102 ^02,112 /'02,32 /'02,92 /'ll2,32 /'l22,102 /'l22,ll2 ^122,32 /'l22,92 /'92,102 



Pn2,i02-/'92,32 Ix\/I«!ei 



+ A 8 + 



2 /'l2i2,4i2 P\2x2Jx2 ^'9i2,12i2 P9i2,4i2 ^'9i2,ei2 /'9i2,/i2 ^'di2,12i2 /'di2,4i2 /'di2,ei2 



-/'Ii2,/i2 - Pei2,4i2 " /'ei2,/i2 ) I >< a/ '^^H^i^'') 



|8 



2 \ ^^45:3i5 ^45-3i5 ^45:Ci5 ^^5,6^5 ^d[5,f[5 ^e[5,c[5 ^f[5,3i5 ^f[5,e[5 ^a[5,3i5 ^a[5,c[5 



2^ ^ _ p2^ ^ 



+ A 8 + 



/'l34,184 Pl34,C4 PW4,n4 /'l64,184 ^164,64 /'64,184 ^64,C4 ^04,134 Pa4,lS4 Pa4, 



Pa4,C4 



|8 + ^ ( ~ /'l3i4,ai4 ~ ^13i4,/ii4 ~ Pl3i4,hA ~ /'/ii4,ai4 ~ Pji4,ai4 ~ Pj^A^hiA ~ Pjii,h^ ~ ^A;i4,13i4 ~ PkiiM^ ~ PkiAJiA 



pUma - PhA,a,A ) I X A/ \det[E^{u) 



+ 



2 2 2 2 2 22222 

+ - ( - Pl75,65 ~ ^175,135 ~ /'l75,r5 ~ ^135,65 ~ ^P5,175 ~ ^P5,65 ~ PP5A'^ ~ Pp5,r5 ~ PqsfiB ~ ^95,135 



PLr3, \ x\ \det{E^{u) 



(3.38) 
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Expanding the square roots (see Section I2.2p and considering just first-order terms we obtain 



Vr = 6\I\ det {E-{u)) I 2 ( 9 + ^ X 1 ( - pl^^^.^ - pl^^^^ - pl^.^^^.^ - p^^^^g^ - pl^^^^^.^ - p\^^^^.^ - pl,.^^^.^ - p\,^ 



2 2 2 2 2 

^'173, 133 ~ ^173,23 ~ Pl33,83 ~ /'133,23 ~ Pl4l3,21i3 



2 2 
/'9l3,14i3 ~ /'3l3,21i3 



^4i3,14i 



2 2 
/'Sia.Qis ~ /'5i3,21i3 



2 2 2 2 

/'5i3,3i3 ~ P^i-iMi-i ~ /'Oa.lOa ~ /'02,ll2 



2 2 2 2 

" /'3l3,9l3 ~ /'4l3,9l3 ~ /'4l3,21l3 ~ /'4l3,3l3 
2 2 2 2 

^02,32 ~ /'02,92 ~ P\\i,Zi ~ Pl22,102 



2 2222 22 2 2 22 

/'l22,ll2 ~ PVilM ~ /'122,92 ~ /'92,102 ~ Pll2,102 ~ /'92,32 ~ /'l2l2,4l2 ~ /'l2i2,/l2 ~ /'9l2,12i2 ~ P9l2,4i2 ~ P9l2,ei2 

/'9i2,/i2 - ^rfi2,12i2 - ^^12,412 - ^rfi2,ei2 " /'rfi2,/i2 " ^ei2,4i2 " ^ei2,/i2 " ^^c'^S.SiS ~ ^45,315 ~ ^3;5,ci5 ~ ^^S.e'.S 



2 _ 2 



2 _ 2 _ 2 

^X5,3i5 ^X5,e'i5 ^45,3i5 



2 2 2 2 2 2 2 

/'l64,184 ~ /'l64,64 ~ /'fe4,184 ~ /'fe4,C4 ~ Pa4,134 ~ /'a4,184 ~ /'a4, 



2 



Pa[5,f[5 ~ ^134,184 ~ Pl34,C4 ~ Pl64,134 

2 2 2 2 

" Pc4,164 ~ /'l3i4,ai4 ~ Pl3i4,/ii4 ~ /'l3i4,/ 



2 2222 2 22222 

Phi4,ai4: ~ Pji4,ai4 ~ Pjj^4,hi4 ~ /'ji4,Zi4 ~ Pfci4,13i4 ~ Pki^M^ ~ Pki4,ji4 ~ Pfci4,h4 ~ P/i4,ai4 ~ PlTsfis ~ ^175,135 



2 2 

Pn5,r5 ~ ^135,65 



2 2 

Pp5,l75 ~ Ppo,' 



2 

Pp5,q5 



f(n-l) 



2 _ 2 _ 2 _ 2 

Pp5,r5 Pqsfib PqsA^s Pq5,r5 



Prs, 135 



det 



5X- 



253^|det(i?»)|(9 + i(-l 5: pL,)) 



det 



6{s, u^,u'^] 



6{s, v} ,v? 

(3.39) 



Here pi^^k^ represent the off-diagonal entries of the matrix ^/~A ^ (see Section lA. 2. 2p which denote the value 



of the term 



for the edges Ci^i and ei^k incident at the vertex Vi. 



The term proportional to p^ ^ represents the higher-order corrections to the expectation value of the 



volume operator. Each term p^^.J^. is proportional to Cx j for /' < / where C is a constant that depends on the 
Euler angles we choose. It follows that only the zeroth-order of the expectation value of the volume operator 
for the 6-valent graph is rotationally invariant, up to embeddings with measure zero in SO{3). However, 
only embeddings which have measure zero in SO (3) (when the edges are aligned to the plaquettes) give the 
correct semiclassical limit. 



3.3.3 Expectation value of the volume operator for a translated 6-valent graph 

In this Section we will calculate the expectation value of the volume operator for a translated 6-valent graph. 
To make the comparison as accurate as possible, we translate the 6-valent graph by a vector with more or 
less the same properties as the vector with respect to which we translated the 4-valent graph, namely: 

1) I > ex > ej/ > > 

2) '-^<e. + e, 

3) ¥ - f < 6. 

4) K\- \yi\ > nil - nil for all |T4^| > \V^\ 
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The coordinates of the translated vertices are 



= {e^,ey,e,) (3.40) 
5V2 



^2 



+ + (3.41) 

= {e.+ '-^-'-,ey-'-,6V2 + e.) (3.42) 

[6V2 + e^,ey-5,5V2 + e,^ (3.43) 

(e,,e^-5,e,) (3.44) 

+ + (3.45) 

Vi', = (5V2 + ^ + e,,e,-^,e, + ^) (3.46) 

^15 = (-^- Y + ^-'^J^-'^'^^ + ^O ^^-^^^ 



14 



Similarly as for the aligned case we have 

\l iff^e, =0 \l iff^e, =0 

F,^ = ^ ^^^'^•^^ Fl = X Ve.G7 (3.49) 

[00 iff^e, =0 ^ ^ 

The expression for the matrix \fA ^ is given in Table IA.2.31 in the Appendix. 
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Given the conditions above we obtain the fohowing value for the expectation value of the volume operator 
as applied to the periodicity cell: 



V; 



R 



'|8 - Ml! _ ^ _ V2i?2 - A',BoB', - - 2{B',Y - BoCo - - ^ - 



V2 



2^2 



— 4 i7! 2 4 i7! ^^l^yl^^H^^-^^^JI 

+ V ' ^ ^-A,B,B, - -^-^ - 3B, - 2V2B, - (C,) - 



2J2 



Cl/2\x ^j\det{^E-{v))\ + ^\8 



3(4)^ 5{B',)^ AsjB',)^ B'.Bs 3A',B',B^ 



s/2 



"^3-^3-^3 od2 o,/^o2 ^3-^3 ^3-^3 ^'s^sC'^ B'^B^C'^ SB'^BsC'^ 3(^)2 ll^^J pa(„\\\ 

, / 3(4)2 A4(^;)2 , , 34^:^4 g,/^p2 ^;^4c; 3g;g4g; 

+ y|8-^^---3(B4) -^^1 4^4^4 - - 354 - 2V2i?, - ^ 



3(c;)2 4c;c4 ci 



2V2 



3BI 13(C;)^ 4(gD^ ^(Cs)^ J?5(0;)^ B5(C;)^ (4)"g5 34c;C5 (C^j^Cs 5C|, 

-D5<-^5 " 



4^2 



2V2 



4^2 



4\/2 



4^2 



8^2 



X sl\del(E^{u))\ + J\8- - - A[,B[,B,s 



^13-^13-^13 

V2 



13 



'13 



4^2 



2V2 



A' R' R 348-^18-^18 qR2 9,/^r2 ^18-^18 C^s 3^i8-^18<^18 3(C^8)2 (^i8)^Ci8 48^^18 '^18 

^^^1 SB,, - 2>/2i5,s ^ j-^= J j-^ 



(^15)^-^15 ^/T^A o 0R2 ,/oh2 45^15^15 -Sig^lsClS -6^5^15(^15 3 ^15^*15 -^15^15 | 



X W|det(i?^(n))| + J|8-^ 



2(s; j2 _ 4^b;45i4 
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x/2 



^2^2^ 



3(c;4)^ 



-Si4(C;^4)2 3i3i4(C;^4)2 A^^B^^Cu ^ A^^C^^Cu B^^Ci^Cu B^^C^^Cu {C-^^'^Cu 3(72^ 

-D14(-^14 — 
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2\f2 



4^2 



8^2 



(3.50) 
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Expanding the square root and considering only first-order contributions (see Section 

det 



we obtain 



VR=6\l\det{E-{u 



6X1 



3(4)^ 



1 



2 9 + — 
16 



5Bi 



y/2Bl 



o; ~ BqCq 
KB's? 



2V2BI 



2\/2S| 
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Al - V2A5B5 - 
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13(^5 



^ - ^{B[,r' 
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B5C5 

18 



3(i?;)^ 

5Cl 



3Bl 



3iA 



13) 



2V2BI 

n{B 



3(Q 



3^0^ 



B,Bs 
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4 



3(^5^ 



13) 



3B 



2V2B^ 



3(c; 



18 



is; 



3^13 



2V2B^ 



3{C 



13) 



13 



2A 



15 




V2B^ 



15 



C 



15 



515C15 



3Kl4)^ 



2(5 



14 



14; 



\/25i\ 



3(c; 



u) 



(3.51) 



where the terms 



A' 



6i^V2 + l) 
{\yv, \ - mil) 



6{yJ V2 + I) 



Bi 



6UV2 + 1) 



i\xv, \ - rijl) 

{\yv^ \ - mji) 



(3.52) 



are the off-diagonal matrix elements of yA (see Section lA.2.3p . The quantities xy. , yy., zy. represent are 
the X, y, ^-coordinates of the vertex Vi respectively. 

In the last line of the equation above we have performed a Taylor expansion and factored out the term 



det (^Ej(u)), since we can assume that, although it is vertex dependent, the values of this term to first- 
order will be the same for each vertex. Due to the appearance of the terms ()3.52p . which are proportional 
to the Euler angles, equation (j3.50p is translational invariant (up to embeddings of measure zero in SO (3)), 
only at zeroth-order. 



3.4 Analysis of the expectation value of the volume operator for 8-valent graphs 

In this Section we will calculate the expectation value of the volume operator for an 8-valent graph. As in 
the case of 4- and 6-valent graphs, we first consider the non-rotated graph. We then analyse the rotational 
and translational dependence of the expectation value by performing a rotation and then a translation of 
the graph; we then repeat the calculation. It transpires that, even for the 8-valent graph, the off-diagonal 
elements of the matrix have non-trivial contributions that cause the expectation value of the volume operator 
to be translationally and rotationally dependent, for higher orders than the zeroth-one. 



3.4.1 Expectation value of the volume operator for a general 8-valent graph 

As in the previous cases, we take the (0, 0, 0) point of the lattice to coincide with the (0, 0, 0) point of 
the plaquette, and each vertex to be symmetric with respect to the axis. The coordinates of the vertices 
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comprising the periodicity cell are the following: 

Vi = 
= 



-25 
0,^,0 



12 



^4 



(3.53) 
(3.54) 
(3.55) 
(3.56) 



For the 8-valent lattice we choose Vi as our reference vertex. The allowed value for its x-coordinate is 
nl < I^Val < ra/ + |, where n = ["^]- Similarly to the cases of 4- and 6-valent graphs, the allowed positions 
of the remaining vertices in the periodicity cell can be computed from the allowed positions of Vi. Because 
of the geometry of the 8-valent graph we obtain 



rj^{x,y} _ p{x,y} _ ^{x,y} _ ^ 



(3.57) 



This results in the following values for the terms 



as computed for the above five vertices. 



V^ 



^9 



Vl2 



all-^ 







14 terms equal to /3; 12 terms equal to 2/3 



8 terms equal to 2/3 



8 terms equal to 4/3; 4 terms equal to 2/3 



Ve 1 term equal to 6/3; 9 terms equal to 2/3; 12 terms equal to /?; 2 terms equal to 4/? 



Here /? := (-^ — n/) -^^ an d it is proportional to the off-diagonal entries of the matrix ^/A whose values 
are given in Section [A.3. II in the Appendix. 

The expectation value of the volume operator is: 



5{s, ,v?) 



|2 (^V^,^\det{E'^{u)^ I + ^\det(^E^{u)'j \ x V|32 - 60/3^ - 28/?3 



+ J\det(E^{u)]\ X V|32 - 16/32| + J \det(E^{u)] \ x y/\32 - 14462 _ 64/3' 



+ ^ |det (S«(n)) I VI 32 - 104/32 - 32/?3| j 
By approximating to first-order (see Section 12. 2p we obtain 

V„=U'Jl,,t(E-iu))l\mi^^^)\ (5+1x1(324/3^) 



(3.58) 



(3.59) 



In this case, the deviation from the classical value of the volume of a region, R, is of the order four, even to 
zeroth-order in 1/5. 



3.4.2 Expectation value of the volume operator for a rotated 8-valent graph 

We will now analyse the expectation value of the volume operator for a rotated 8-valent graph. In order 
to make the comparison with the 4- and 6-valent graphs as accurate as possible, we will rotate the 8-valent 
graph by the same amount the other valence graphs were rotated. It follows that the angles of the edges 
incident at Vq will satisfy the following conditions: 
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1) < ^0,2,^eo,8 < f > ^ < %,^0,7 < 27r f < 6/0,1,^0,6 < VT and VT < ^0,4,^0,5 < X 

2) ^ < (t>eo.2^4'eo,3 < 27r-sin"^i, 7r + sm~i(i) < cl)eoj,(j)eo,8 < X'T" < '^eo,6 > (^ecs < vr-sm-^(i) and 



sin < 0eo_4,(/)eo_^ < f 

The angles for the co-hnear edges are defined through the formula 9p = 180° — 9p „. and (t>e = 

o o o c ^colhncar " ^ 

180° - ^ecoiiinear fespectively. It follows that the edges eo,i, eo,2, eo,3, eo,4, £9,5, eo,6, eoj and eo,8 he in the 
octants B, A, D, C, G, F, H and E respectively. The coordinates of the rotated vertices are 

Vl = ((--Rll--Rl2 + fll3)-^,(--R21--R22 + fl23)-^.(-B31--R32 + fl33)-^) (3.60) 

Vi = (-2if,2-|,_2ii22-|,-2fl32^) (3.61) 

3 5 S 

VI2 = ((-i?12 + 2i?i3)^,(-i?22 + 2i223)^,(-2i?32 + 2i?33)^) (3.62) 



V! 



({—Rii — Ri2 + '^'^i3)"y^i (— -R21 — -R22 + '^'^^s)-^) (— -R31 — -R32 + '^'^33-''^) (3.63) 



As was done for the 4- and 6-valent graphs, in order to carry out the calculations for the expectation value 
of the volume operator we would have to specify a particular combinations of angles satisfying conditions 1) 
and 2) above. However, all combinations satisfying 1) and 2) above lead to the same value, in zeroth-order 
in |, of the expectation value of the volume operator. Rotational dependence will only appear for higher 
orders in |. Moreover any sub-case of 1) and 2) will lead to the same couples of edges commonly intersecting 
surfaces ^ in a given stack. Therefore, to leave the result as general as possible, we will not specify a 
particular sub-case of 1) and 2) but simply derive a general expression for the expectation value of the 
volume operator given conditions l)and 2). 

The expectation value of the volume operator for a periodicity region, R, is then computed as 



1 / 1 



det (£«(«))! 4(^5 + -(-- ^ al)i.t[j^^)'^ (3.64) 

where the terms Uij are the off-diagonal entries of the matrix ^/A ^ (See Section IA.3.21 in the Appendix) 
Evidently, the higher-order corrections are angle dependent, while the zeroth-ones are not. Therefore, 
as for the 4- and 6-valent case, the expectation value of the volume operator for the 8-valent graph is 
rotational invariant, in zeroth-order up to measure zero in 50(3). However, it does not reproduce the 
correct semiclassical limit. 

3.4.3 Expectation value of the volume operator for a translated 8-valent graph 

We now consider the translated 8-valent graph. As for the 4- and 6-valent graphs we choose the following 
conditions on the components of the translation vector: 

1) 6 > > ej, > > 

2) iVfci - \VI\ > nil - nil for ah \Vl\ > 
Similarly as for the aligned 8-valent graph we have 

rp{x,y} ^ p{x,y} ^ 2{x,y} ^ ^ Ve^ G 7 (3.65) 



50 



The coordinates of the translated vertices are 



12 



V 



{(-X 





(3.66) 


5 5 5 \ 


(3.67) 




(3.68) 




(3.69) 


5 5 5 \ 
73^'"'"73^''' 73^'V 


(3.70) 



The expression for the matrix ^/~A ^ is given in Section IA.3.31 in the Appendix. 
The value obtained for the volume of a region R is: 



Vr 



Is' 



det 



6X§ 



6{s, tt^, li^ 



|4(32 - 4^2 _ 12^2 _ 2AoBl 



2AoC^ - QBoC^ 



BoDoEo _ ^2^^ _ clF, - F2 ) | x ^/ \det ( i^«(n) 



+ W |4 32 - l&Al - 2A\ - 24Bl - SAiBf - BfDi - Df - AjEi - BfEi - Ef - SBfFv 



AiDiFi 



3F2-M? 



X J\det[E''{u) 



+ y |4(32 - 2451 - ABlC - 8C| - QBID, _ ^ _ ^ _ - ^ - ^ - BjE,- 



-^^-Fl)\.Mdet{E^{u) 



+ a/|4 32 - 24^22 - 6^22^12 - 12^22 - 2^2^Ci2 - 2522C12 - 4C22 - Al^Di2 - Bl^Di2 - Dl^ - Al^Eu - Ef. 



12 



-3^2^F,2 - _ - ^) I X ^\det{E^i^ 



, 2M2 43 R2 

+ i/|4 32 



6AeS2 - 2AlCe - 2Cl - BID, - Dl - 2AlE, - BlE, 



3E^ 



' 3BlFe- 



3Fi-^]\^^\\det[E<l{u) 



(3.71) 
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Considering only contributions of first-order in | (see Section 12. 2|) we obtain 



VR = 46\\det(E-{u) 



det 



6X- 



1 + J_f _4^2_i2R2_1^ 

32 V 2 



3D^, - El - 



+ A 1 + 



32 



IQAl - 24^2 -Dl- -El - 3Fl ) + \ 1 + 



1 

32 



+ A 1 



32 



24^2^ - 12i?22 - 4C22 - Dl, - El, 



1 

32 



A5-'\j\det(E<^{u] 



deti 



5X% 



X |5 + ^x^( -4yl^-12i?^ 



45C72 



3D?, - El - F^ 



IQAl - 24BI -DI- -EI - 3Ff 



2ABi 



M 3El 



- 24^?2 - 12S?2 - "^Cl, - Dl, - El, - 3Fl, 



29A? 43^2 



2^2 - Di 



3^2 



3f: 



(3.72) 



where Ai, Bi, Ci Di, Ei and Fi are the matrix elements of \/]4 ^ as defined in Section lA.3.31 of the 
Appendix. Again, translational invariance holds only at zeroth-order up to measure zero in SO{3). However, 
it does not reproduce the correct semiclassical limit. 



4 Summary and Conclusions 

We have shown that if we use semiclassical states derived from the area complexifier, then we do not obtain 
the correct semiclassical value of the volume, operator unless we perform an artificial re-scaling of the 
coherent state label and we restrict our calculation to the following special cases: 

1) The edges of the graph are aligned with the orientation of the plaquettes (6-valent graph). 

2) Two or more edges lie in a given plaquette (4-valent graph). 

3) One edge is aligned with a given plaquette while a second edge lies in a given plaque (4-valent graph). 

However, such combination of edges have measure zero in S0{3). For embeddings whose measure in S0{3) 
is non-trivial we do not obtain the correct semiclassical behavior for the volume operator for any valence of 
the graph. 

This result suggests strongly that the area complexifier coherent states are not the correct states with 
which to analyse semiclassical properties in LQG. Moreover, as previously mentioned, if embedding inde- 
pendence (staircase problem) is to be eliminated, area complexifier coherent states should be ruled out as 
semiclassical states altogether. 
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A Edge metric calculations 



A.l Edge metric components for the 4-valent graph 

In this Section we list the values for each of the terms ''"^ which comprise the matrix ^TA ^ that appears 



in the computation of the expectation value of the volume operator. The method for computing such terms 
was described in Section [3.11 with the aid of a two-dimensional example. It is worth recalling that the term 
te^ej = ifie, + *eie, + ig.e. represents the total number of plaquettes in each direction that are intersected by 



are the total number of plaquettes in each direction intersected by the edges and e,- respectively. 



both the edges Cj and Cj. On the other hand, the terms te^ = tg,, +tei +tl. and te, = t^- +tej +^ej represent 



A. 1.1 Rotated 4-vaIent graph 

In this Section we will list the values of all the terms, ^=i=, that appear in the matrix \/A ^ for the 

rotated 4-valent graph. This is a 15 x 15 matrix labelled by the edges belonging to four vertices, as for the 
aligned case. At this point, we recall the meaning of the following symbols 

cost^p. smUp. 
F^ = |cot0e | py = \ tan Oe ] 

Ci \ Cj I Ei \ till 

Tf. = I tan (pei cos 6*6 . | T^^ = \ tan sin | 

*ei = tf. + t^. +tl.= 5 sin(/)e. cos6'ei + 5 sin(/)e. sin6'e. + 6 coscpei = asiS (A.l) 
which represent the angular dependence of the terms ^t^, ^t^, ^t^, ^t^ , ^t^ and ^t^ respectively (see Section 

EH). 

The arrangement of angles for which we choose to perform our calculations will be such that the edges 
incident at vertex Vq — eo,i, eo,2) eo,3 and eo,4 — lie in the octants F, A, H and C respectively. However, as 
mentioned in previous Sections, in order to compute the specific value for the terms t^.^ej, we would need 
to choose a particular value for each of the edge's angles. Nonetheless, we will refrain from doing so here, 
since we would like to leave our result as general as possible. What we will do, instead, is to write the terms 
tei^ej as functions depending on the angles of the edges and on the position of the vertex at which they are 
commonly incident. Specifically, we recall from Section [3.11 that the number of plaquettes a given edge 
intersects in direction z is given by 

tl:=^tlnyti:=Ti (A.2) 



where 



and 



(xvj — nxy.l)Z^. iff the edge points in the negative x-direction 
(nxv^ + ^ ~ )Zg, iff the edge points in the positive x-direction 



{{xvj — nxY.l)Zi^ iff the edge points in the negative z-direction 
[rixy l + / — xy.)Z^^ iff the edge points in the positive z-direction 

Similarly, for intersections in the x- and y-direction we have the following: 

tl^ r\y tl := Tl (A.3) 



where 



yti 
til 



(uVj — Uyy l)F^, iff the edge points in the negative y-direction 
{uyy l + / — yy.)F^. iff the edge points in the positive y-direction 
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and 



where 




n. 



iff the edge points in the negative z-direction 
iff the edge points in the positive z-direction 



n. 



l)Fl 



iff the edge points in the negative y-direction 
{uxy l + I — xy. )Fe- iff the edge points in the positive y-direction 



(A.4) 



n 



2V, 



l)Tl 



iff the edge points in the negative z-direction 



[n 



zv. 



I + 1 — zy, )T'e- iff the edge points in the positive z- direction 



From the above formulae it is straightforward to deduce that the precise value for the terms ^t^. will differ 
according to which angles we choose for each edge. However, which couple of edges commonly intersect 
surfaces in a given direction only depends on the octant in which the edges lie. Since we do not want to 
restrict our calculations to any specific angle, we will simply define the terms tei,e as follows 



Tl^ n Tl^ + Tl n Tl + Tl^ n vi^ -.= E,^ n 



and 



n: 



(A.5) 



(A.6) 



To simplify the notation, we denote the edge, e^, going from vertex Vt to vertex Vj as e^j. We can then 
write the terms above in the following simplified way: a^^. ,,6^ ^ = ctki,ii- 



Given the above, the non-zero entries of the upper-half entries of v A for the rotated 4-valent graph 



-1 



are listed below. 










(^)5,2 = -^«72,02 


(^^)6,2 = - 


2 '-'^62,02 


= - 




(^)6,5 = -^"62,72 


(VI)7-,5 = - 


"2*^52,72 




2*^52,62 


(^)9,8 = -I"58,128 


(^^)ro!8 = 


-^ail8,128 


{^)ll]8 = - 


-|ai08,i28 


(^^)k),9 = -I"ll8,58 


i^)n]9 = 


-^"108,58 


(^)n!io = 


— |ai08,ii8 


(^)l3,12 = -I«15i3,3i3 


(^)m!i2 = 


— 2CKl8i3,3i3 


(^^)r5!i2 = 


"2*^1913)313 


{'/^)lll3 = -^"I8i3,15i3 


(^)r5!l3 = 


-- -^ai9i3,i5i3 




-|ai9i3,18i3 


Here (\/]4)j j indicates the i-th column and the j-th 


row entry of the matrix 







A. 1.2 Translated 4-valent graph 

In this Section we give the matrix entries of ^/~A ^ obtained after translating the 4-valent graph. As for the 
aligned and rotated graphs discussed above, there are four vertices in the periodicity cell. The values for 
the terms te^ej = t^^ej +*eiej +^eiej ^s obtained for each vertex are given below. The (Gauss brackets) terms 



"TJ' 



and Pi 



represent the number of stacks which each edge Cj intersects in the x-, 



y- and z-direction respectively (see Section [3T]) . 



1) Vg" = (ex, Cy, £2) with ex> ey> 



2) K 









*eo,ieo,3 ^ ^60,260,4 ^ (^S/ ~ "'■oO 


*eo,3eo,4 ~ *eo,ieo,2 ~ (^^ ~ Pol) 


~ ^£0,360,2 ~ ^eo,ieo,4 ~ (^2 ~ PoO 



V3' 
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62, 662, 5 


= *62,062,r = (l^^-:^|-^20 


*e2,oe2,B ~ *e2,7e2,6 ^ + - P2I) 


^62, 062,6 ~ ^62,662,7 ~ (l^a; !^ 1 









68,1068,11 


— f 2 — _|_ _ 

~ ''68,568,12 ~ ^ 




4.x 4.x 

68,568,10 68,126 




''68,1168,5 ''68,1068,12 





5)^i'3 = (6.,e, + 2ife,6,-2i^) 









^613,3613,15 ~ ^613,19613,18 ~ ~ ^IsO 


fX _ fX _ n 

613,19613,15 613,18613,3 Vl^^: 


-2^|-P130 


fV - fV _ 

''613,15613,18 ''613,19613,3 

= {f-x - ni^l) 



The upper-half entries of the matrix \/~A are hsted in the fohowing table. It should be noted that the 
lower indices refer to the matrix entries, while the upper ones indicate the edges that we are considering. 
Therefore (V^ )2 V^" denotes the matrix entry at row 1 column 2, which is the term = correspond- 

ing to the edges eo,i and eo,2, which join vertex Vq to Vi and Vq to V2 respectively. 



-1n1o,2o 

>2,1 
"^\2o,3o 
^3,2 

-Ixl9l3,18l3 _ 1 A 

79,8 — 2^ 

"I\18l3,15l3 \ i 

jl0,9 --2^13 



(VI" 
(VI" 



13 



^108,118 

n3,i2 

^118,128 

n4,i3 



2^ 
i^8 



where Ai = - riil)^, Bi = (|?/y, 

In all the above formulae, the term 



(VI" 


■-'^\loi3o _ 
^3,1 — 




(VI" 


■l\lo,4o _ 
>4,1 — 


— ^Cq 


(VI" 


-l\2o,4(, 
M,2 — 




(VI" 


-1\02,72 
l5,2 — 




(VI" 




\M 


(VI" 


-l\3o,4o 
M,3 ~ 




(VI" 


-^72,62 

— ~ 




(VI" 


-1x52,62 

)7fi — 


— \B2 


(VI" 


-Ia9i3,15i3 
J 10,8 


= -^Ci3 


(VI" 


-Ixl9l3,3l3 

^11,8 


= 


(VI" 


"I\18l3,3l3 

yii,9 — 


-5C13 


(VI" 


-I\15l3,3l3 

711,10 


= -I An 


(VI" 


-1x108,128 _ 

714,12 — 




(VI" 


-1x108,58 _ 

715,12 - 


- —\c% 


(VI" 


-1x118,58 _ 

n5,13 ~ " 


-5C8 


(VI" 


-1x128,58 _ 

715,14 - 





(5eV3 



comes from the fact that te - 



5e V3 for all edges 



A. 2 Edge metric components for the 6-valent graph 

In this Section we give the values of the terms ^J''^"''^ that represent the entries of the matrix VI ^ for 
the 6-valent graph. We do this for all non-rotated, rotated and translated cases respectively. 



A. 2.1 Aligned 6-valent graph 

For the 6-valent graph, the periodicity cell contains nine vertices. Thus the matrix VI ^ is a 46 x 46 matrix 



We will now give the values for the terms 



''k,j 



which we denote as ^ikJk- "^i^^ 
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following short-hand notations 



V2 



— nl 



V2 



— nl 



a 



1 "; 

2 -n / 



6{l + 4 



) 



a 



: -- 



SVV2 + 1 
1 "l 



,5(1 + 



A 



k ■- 



6VV^+ 1 
6VV2 + 1 



1) Vo all terms Tj^^j^^^ = 

2) ^2 



3) ^3 



4) V4 



5) 

6) Fl3 



7) Fl8 



^02,112 
^02, 92 = 

^122, 32 
^32, 92 = 



a 
2 

a 

'A 

_ a 
V2 
aV2 
4 



J- 02,32 = 
^112, 122 

ri22,i02 



a 
' 2 



a 

a V2 



J- 02,102 - 
^112,102 
ri22,92 = 



a 
' 2 

4 

a'V2 
4 



^23,173 = 

^23,163 = 
^173, 153 
^193, 153 



2 



ra4,184 - 
-^04, 64 — 

ri84,i64 

^64,04 = 



-^(72 + 1) 

a'y2 



r23,i33 - -f^i^ + i) 



2 



J- 5l3,9l3 - 

r5i3,i4i3 

^413,2113 
^913,1413 



^418,918 

r4i8,ei8 

■^618,918 
■'^Cil8,/l8 



a 2 

V2 



r5,175 



a 
' 2 

/ 

a 
2 



V2 



^173,133 
^133,153 

1^193,163 

ra4,134 ~ 

ri84,i34 

1^134, C4 ~ 

ri64,C4 ~ 



+ 1) I 



lis 
2 

A3 

2 

/ 

V2 



-a 

_^ 

V2 

/ 

-a 



23,193 



173,193 
'163,133 
'163,153 



r — 

ri84,64 = 
ri34,i64 
rfe4,i64 = 



As 
2 

2 

a'V2 
2 

V2 



2 



~V2 



65,95 



aV2 
4 



-'^5l3,3l3 - 
^413,913 - 
^413,1413 
^313,2113 



_ Q 

2 

_V2a_ 

4 

a'V2 
4 

^ 



2 

V20 



^18 
2 



2 



4is,12i8 

ri2i8,9i8 
ri2i8,/i8 

■'^618, /18 



V2 + 2' 



2 

aV2 
4 



^513,2113 = 
^413,313 = 
^313,913 = 

r2ii3,i4i3 

r4i8,di8 = ■ 
r/is,9i8 = - 

rdi8,12i8 — 
■■- ei8,di8 



a 
2 

2£ 

aV2 
4 

4 



dm 

"V2 

_^ 
/2 
a'V2 



(A.7) 
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15 



9) ^14 



13i4,ai4 



r, y = 

•J15,"i5 


■ 

2 


r, ' = 

3l5,Ci5 


3a'\/2 
4 


''^315, /is 


r, ' = 

.-il5,aig 




r ' = 




r / / = 

•^15 '"^15 


r / f' — 


3a' 


r / , = 




Tc' d' = 




2 




2 


'-15'"lB 






r / / = 


3a' V2 
4 


r / f' — 



3/3 



r/ii4,ai4 

■'^J 14,^14 ■ 
r-1 



^14 
- V2 



14 
2 



' V2 
2 



2 



13l4,/ll4 



r'/i4,fei4 



Ai4 
2 



'^15 
V2 



3£ 
V2 



The upper-diagonal entries of the matrix s/A are hsted below. Similarly, as for the 4-valent case, 
lower indices represent the matrix entries while the top indices indicate the edges we are considering. 
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^fA 

Va 

VA 
VA 

Va 
Va 
v^" 
Va 
Va 

VA 

Va 
Va 

Va 

VA 

Va 
Va 
Va 



-1 



-1 



~i 



~i 



-1 



-1 



~i 



02,112 _ 

7,2 — 

02,92 _ 

11,2 — 

122,32 _ 

9,8 ~ 

32,92 _ 
11,9 ~ " 

23,193 _ 
14,9 ~ 

173,193 _ 

14.12 - 

133,163 _ 

16.13 ~ 

3l3,9l3 

20.13 ~ 

193,163 _ 

16.14 - 

3l5,Ci5 _ 

39.15 ~ 

a4,134 

19.17 ~ 
184,134 _ 

19.18 ~ 

4l8,9l8 _ 
32,18 ~ 

4i8.ei8 _ 
36,18 ~ 

4l3,14l3 _ 
29,19 

C4,fe4 

21,20 ~ " 

»'15,17l5 
23,24 

5l3,21l3 

31,25 

14l3,21l3 

31,29 

9l8,12l8 

33,32 

12l8,/l8 _ 

35,34 

/i8:ei8 _ 
37,35 — 

40,38 — 
d' c 

"15 '"-15 

41,40 ~ 

r' ' 

43,42 — 

ai4,/ii4 _ 
47,44 - 



V2 

/ 

_ ^ 

as/2 
4 

_^ 
2 

2 

a V2 
2 

aV2 
2 
/ 

V2 
4 

-a 
_^ 

V2 



n/2 

4 

ay/2 



P 



2a 
' V2 



a 
' 2 

ay/2 
4 

2 

ay/2 
4 



2 



V2 + 2^ 



^15 

f 

■ V2 



02,32 _ 

9,2 


a 
2 




122,122 _ 

8 7 — 


r 

a 




122,102 _ 
iU,o 


a y/2 
4 




23,173 _ . 

12,9 


2 




23,163 

16,9 ~ " 


p' 
2 




173,153 _ 
15,12 


-/?' 




3l3i4l3 

17,13 — 






3l3 ^21 13 

21,13 ~ 






153,163 _ 


a; 




17,15 


V2 


3l5>/i5 _ 
40,15 


-Aim. - dp 

V2 ' 




a4,C4 


ay/2 




20,17 ~ ~ 


4 


184,64 
21,18 ~ " 


-b' 


\, V ^ 


4l8,12l8 

33,18 ~ 






134 ,C4 

20,19 ~ " 




fx/I'' 


4l3,9l3 


a y/2 


fx/I"' 
I, V ^ 


30,19 ~ 


4 


C4,164 

22,20 


a 
2 


(x/I"' 


5l3,14l3 
29,25 ~ 


{^^A~^ 


95,65 _ 
28 27 — 


-/3(\/2 + l) 


(x/I-' 


13i4,ai4 _ 
43 29 ~ 


2p 
V2 




9l8,/l8 

34 32 — 


p 

9 




12l8,lil8 
36,34 ~ 


/o 
v2 




(^18,618 


a y/2 


(VI-' 


37,36 — 


4 


' ' / 

ei5>/l5 _ 


3o' 


(VI-' 


41,39 ~ 


2 


'^'iS'/w _ 
42,40 ~ 




(VI-' 


fci4,ii4 


^14 


(VI-' 


45,43 ~ 


2 


'l4,jl4 
46,45 ~ 







02,102 _ 

10,2 — 
112,102 _ 

10.7 ~ 

122,92 _ 

11.8 ~ 

23,133 _ 

13.9 ~ 

173,133 _ 

13.12 — 

153^-1 
715,13 

3l3>5l3 _ 

18.13 ~ 

193,153 _ 

15.14 — 

315,^15 _ 

38.15 ~ 

04,184 

18,17 ~ 
(14,64 

21.17 ~ " 
184,164 _ 

22.18 ~ 

4l8,dl8 _ 

35.18 ~ 
134,164 

22.19 ~ 
4l3,21l3 _ 

31,19 

64,164 

22,21 ~ ■ 

5l3,9l3 _ 
30,25 ~ 

14l3,9l3 _ 

30,29 

13l4,/ll4 

46,29 

9i8,ei8 _ 
36,33 ~ 

/l8,rfl8 _ 

36,35 ~ 

'^15''^15 

39.38 — 

'^15''*15 

42.39 — 

'^15''^15 

43,41 — 

fcl4,il4 _ 
46,43 ~ 



a 
' 2 
ay/2 
4 



a 



V2 



2 



2 



a 
' 2 



Al5 

2 

_ 2£ _ 

a y/2 
2 

ay/2 
2 

"^18 
V2 



2§_ 
■V-2 
— 

2 

V2 ' 
Ai4 
y/2 



y/2a 
4 

M 
2 



2 



-/?(V2 + 1) 



Al5 

2 

Za y/2 



^14 
2 



A. 2. 2 Rotated 6-valent graph 

In this Section we give the values of ah the terms 



V 



that appear in the matrix \/ A for the rotated 



-1 



6-valent graph. In order to carry out the calculations we will rotate the graph such that the edges incident 
at vertex Vq: eo,i, eo,2) eo,6) eo,7) eo,8 and 69,17 lie in the octants ff, -D, G, F, A and B respectively. As 

for the aligned graph, we consider the periodicity cell composed of nine vertices. Therefore VI will be a 
46 X 46 matrix. ^ 

In what follows the terms Pi^A^, represent the values for = and are defined as for the 4-valent 



case (see lA.l.T]) 

Given the above, the upper-half elements of the matrix 



A 



V*^M*^*^7 



are: 
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V 


-1 

20,2 " 




— — r)r\ n 
2K02,92 


(y/A: 


-1 

22,2 " 




— — /)n 1 1 

2^02,112 


(Va' 


-1 

24,2 " 




2^02,32 


V 


-1 

25,2 " 






v ^ > 


-1 

8,7 ~ 




— 1 c 

2 Pa4,lo4 


(Va' 


-1 

9,7 ~ 




2Pa4,134 


V 


-1 

10,7 " 




2Ka4,C4 


{y/A: 


-1 

11,7 " 




2/^04,04 




-1 

9,8 ~ 




2^134,184 


V 


-1 

11^8 " 
35,8 " 




2^04,184 


v * > 


-1 

12,8 " 




~ 2/5164,184 




-1 

21,8 " 




— — Dn A 
2F9i8,4i8 


Va' 




— — Ha a 

2Ka!i8,4i8 




-1 

36,8 




— — Da 

2F4i8,ei8 




-1 
10,9 




— — 1 Q 

2FC4,134 


V -'^^ 


-1 

12,9 " 




~ 2/'l64,134 




-1 

19,9 " 




2P4i3,9i3 




-1 

27,9 " 




— — /In A 

2F3i3,4i3 


V 


-1 
31,9 




1 „ 

2/5413, 2I13 




-1 
32,9 




1 „ 

2/^4i3,14i3 


V ^ > 


-1 
33,9 




1 n 

2/^4i8,12i8 


V 


-1 
11,10 




2f04,C4 


{y/A: 


-1 
12,10 




— — 1 a 
2FC4,lb4 


V ^ > 


-1 
12,11 




— — Oi, 1 
2Fo4,184 


V 


-1 
14,13 




— — 1 -7 

2fr5,175 




-1 
15,13 




2f?-5,135 




-1 
16,13 




2 Pr5 ,P5 


Va' 


-1 
17,13 




2/^'^5,g5 




-1 
15,14 




"2/5135,175 


(Va" 


-1 
16,14 




2FP5,175 


Va' 


-1 
18,14 




"2/^65,175 




-1 
17,15 




2/^95,135 




-1 
18,15 




"2/^65,135 


V 


-1 
19,15 




J- /1q c- 

2/^9i3,55 




-1 
27,15 




2A^3i3,5i3 




-1 
31,15 




2F21i3,5l3 


V 


-1 
32,15 




~ 2/'l4i3,5l3 


i\/A' 


-1 
17,16 




~ 2/'95,P5 


v ^ > 


-1 
18,16 




— — Da ^ 
2Fb5,p5 


V 


-1 
18,17 




2^65,95 


v ^ > 


-1 
27,19 




— /In n 

2/^3i3,9i3 




-1 
32,19 




~ 2/5l4l3,9i3 




-1 
26,20 




1 „ 

2f 122,92 




-1 
27,20 




2^32,92 


(Va' 


-1 
28,20 




~ 2/'l02,92 


Va' 


-1 
33,21 




~ 2/^1218, 9i8 




-1 
34,21 




2/^/18,9l8 




-1 
36,21 




2Fei8,9i8 


V 


-1 
23,22 




"2/^122,112 




-1 
24,22 




2^32,112 


(a/A' 


-1 
25,22 




2/'l''2,ll2 


V 


-1 
24,23 




1 „ 

2/^32,122 


i\/~A' 


-1 
25,23 




2/'l''2,122 




-1 
26,24 




~ 2/'l'''3,23 




-1 




Z ' ^'^o ,^3 


{Va: 


-1 




— ^PW-i 2-x 
2r ^-^S ,^3 




-1 

Qn OA 




— ipi6^i 2i 
2 1 J^^3 ,^3 


VA] 


-1 

27 26 




-2^133,173 


{VA] 


-1 

28 26 




-2^193,173 


{Va: 


-1 
29 26 




-2/5153,173 


Va] 


-i 

29,27 


= 


-5^153,133 


{VA] 


-i 

30,27 





-5^163,133 


{Va' 


-1 
31,27 





— 5P21i3,3l3 


Va] 


-1 
29 28 


= 


— ^Pl53,193 


{Va] 


-1 
30 28 


= 


— ^Pl63,193 


{Va: 


-1 
30 29 


= 


— |P163,153 


^: 


-i 

38,29 


= 


1 „ 

2/^d' ,3i^ 

15 ' -'-'^ 


{VA] 


-i 

39,29 


= 


1 n 

2 Pc, c ,3i 5 

1-15, '-'ID 


{Va: 


-1 
40,29 


= 


1 n 

2^X5.315 




-1 
41,29 


= 


_i , 

2 "a-^g,3i5 


{VA] 


-1 
32,31 


= 


-2Pl4i3,21i3 


{Va: 


-1 
42,32 


= 


— 2/''^14,13l4 


VA] 


-1 
43.32 


= 


~ 2/''ii4,13i4 


{VA] 


-1 
45,32 


= 


~2Pkl4,13l4, 


{Va: 


-1 
46,32 


= 


~ 2/''l4,13l4 


Va' 


-1 
34,33 


= 


~ 2/^/18, 12l8 


{VA] 


-1 
35,33 


= 


~ 2/''^l8,12i8 


{Va: 


-1 
35,34 


= 


~2P'^wJw 


^: 


-1 
36,34 


= 


~ 2/'ei8,/i8 


{VA] 


-1 
36,35 


= 


1 n 

2/5ei8,di8 


{Va: 


-1 
38,37 


= 


~\p<,v, 

2^a'i5.e'i5 


Va] 


-1 
39,37 


= 




{VA] 


-1 
40,37 


= 


"f^/(5.45 


{Va: 


-1 
41,37 


= 


Va] 


-1 
39,38 






{VA] 


-1 
40,38 




-^/'/(5.<5 


{Va: 


-1 
41,39 






^: 


-1 
41,40 




1 


{Va: 


-1 
45,42 




1 n 

2/5/ll4,ai4 


{Va: 


-1 
46,42 




~ 2^il5>ai5 


VA] 


-1 
47,42 




{Va] 


-1 
44,43 




1 

2/5jl4,/ll4 


{Va: 


-1 
45,43 




1 n 

2Pji5,ki5 


^: 


-1 
45,44 




~2PkiA,ju 


{Va: 


-1 
46,44 




~2Ph4,jl4, 


{Va: 


-1 
46,45 




~2PhbM5 



A. 2. 3 Translated 6-valent graph 

We will now give the values for the terms tgiej = t^^e + *e,ej + if^e computed for each of the nine vertices 
comprising the periodicity cell for the 6-valent graph. The (Gauss brackets) terms Ui = [-|-], 

and Pi = [-|-] represent the number of stacks which each edge intersects in the x-, y- and z-direction 
respectively (see Section [XT]) . 

1) Vq = (e^, ey, e^) with > ey > 
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„ = (e„ — mnl)Z}! 

60,660,7 ^ y ""yJ"/ fj 

*eo,ieo,7 = - "00^7 
4,260,8 = - "00^8 

4,260,17 = i^y - mol)Zf^ 
4,860,17 = i^y - mol)Zf^ 


4.160.6 = (^^-^'oOrf 

60,260,6 ^'U^y 5 
4,260,1 =(e.-PoOrf 

4.860.7 = ~ P00^7 

4,860,17 = (^y - "ioO^f? 

4,760,17 = (^^-f^oOrf 


Lf^n f>i\ IT — 1 1 Ur\ L ] J- C 



2)F; = (^ + e,,6„^ + 62) 









4,062,10 = (e- + ¥^- "20^0 

4,962,0 = (^2/ - ^20^9 
4,962,10 = - "^20^9 
4,1162,12 = (e- + ¥^-^20^fl 
4,362,12 = - "^20^3 

4,1162,3 = (e. + ^ - nsO^f 


4,1262,10 =(¥^+e^-P20rfo 

^62,962,12 = (^2/ - ?«2i?)-p9 

^62,962,10 = - rn2P)if 
4,1162,0 = (¥^+^--f 20^0^ 
4,1162,3 = (ej/ - "i20^3^ 

^62.362.,, = (ej/ - vn2V)Fi 


^62, 1062, 11 — ( 2 P20^10 
''62,962,3 — \ 2 ' '^Z y2l')-'-g 



3) < = (e. + 56(f - 6, - 1, 2^ + 62) 



(^3> 



— Cf: -I- A (y^ - i 



63,263,13 

63,1763,2 

63.1763,13 

63.I663.15 

63.I663.19 



63,1963.15 



1-^3/1)^1 



n30^f7 
n3l)Zf, 



{\ey-^-msl\)Zy, 
{\ey-^-msl\)Zy, 



63,1663.2 
63,1663,13 



63,1363,2 



63,1563.17 
^X 

63,1563,19 

^x 

63,1963,17 



(|€, + |-m3/|)F| 

= (£2 + 2^-^3/^15 

(k, + |-m3/|)F| 
= (e2 + 2^-p30rf7 

= (|e,-f -mg/DFfg 



''63,1963,2 " 
''63,1563,13 



{e^ + 5,{^-\)-n^l)Ff, 



4 — (2 + Ca;, e^^ 2-2 , ^ 2 









4,ae4,18 = - 41 - mil)Zlg 

4,1364,. = + 2^ - m40^f3 

4,1364,18 = (1^2/ - 41 -m4/)^f8 

4.664.16 = (l^?/ - 41 -mil)Zl^ 

4.C64.16 = (1%- 41 -"140^16 

4.c64., = (e- + 2^-n40^? 


4.a64.18 = (1% -41 -m40i^f8 

4..64.. = (e. + 2^ - p4/)Tf 
4,1864,6 = -41 -^40-^^11^ 
4.1664.13 = (ky -41 -mil)FfQ 
4.1664.C = - - mil)FfQ 

4.1364,. =(e^ + 2^ -P40?^f3 


^64,1664.18 — (^2 +2 2 P40^18 
4,.64,c = (e.+ 2^-^40?;" 







/■^ — 

65,965,6 




- 5e\ - m^l)Zl 


fX 

65,965,6 


= (62- 




= 


(Cx - 




fX 

65,q65,13 


= iez- 




— 

66,665,p 


(Ca; - 


n^i)z; 


^X 

''65.665.13 


= ie,- 


-P5l)Ti 


65,1365,17 


= i^x 


- n^VjZf-j 


fX 

e5,p65,17 


= {ez- 


-P5l)T^ 


65,r65,l7 


-- (ex - 


- n5l)Z^ 


fX 
65,7-66,17 


= {\ey 


- 41 -m^l)F^ 


65,7-65,13 


-- i^x - 


- n5l)Z^ 




= ie.- 


P5IK 



fy 

''65,1765.6 
fy 

''65,re5,9 " 



(ez 



-P5l)Ty 

p^m 
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6) V^^ = 



{\ey -6e\- mi^VjZl 



= 

613,14613,9 

*ei3,5ei3,i4 = (e-^lsO^f 



613,5613,9 



613,3613 4 



613,3613 21 



613 4613 21 



{\ey - 5e\ - mi3/)Z| 
{\€y - Se\ - mizl)Zl 
-- {\ey - 5e\ - mi^VfZl 



4.x 

613,9613,14 



613,4613,14 
fX 

613,4613,9 
fX 

613,3613,5 

613,3613,21 
^X 

613,5613,21 



{\ey-5e\-mizl)F§ 

(e. + ^-pi30rf4 
{\ey - 5e\ - mi3l)F§ 
{\ey - 5e\ - mi3l)F^ 

i\ey-5e\ -mi3l)F^ 

(e. + ^-pi30?5 



^613,3613,9 — {^z + ''2 ^ Pl3i)Fg 
^613,14613,21 = i^Z + ^"2^ Pl3J')Ti4^ 



7) ^1'; = {ex + ^6(2^ + 1), ej/ - I, + ^) 



SeV2 

2 

<5eV2 _ 



6l8,/6l8,e 



618,6618,9 



(|e,-||-mi8/)^| 

{e^ + 5e{2^ + \-nisl)ZI 



6l8,/6l8,9 ^\ y 2 



mi8l)Zj 



618,12618,. (|e^-f |-mi8Z)Zf2 
= (k,-||-mi80^f2 
= (e^ + ,5e(2^ + i - ni8l)Zf 



618,12618,4 
6l8,d6l8,4 



fX 

6l8,/6l8,e 



6l8,dei8,e 

6l8,/6l8,d 
fX 

618,12618,9 
618,12618,4 
618,9618,4 





5e| 

2 1 


- mi8/)F/ 


(e^ + 


2 






5e| 

2 1 


- misl)Ff 




5e| 
2 1 


- mi8l)F^2 




5e 1 

2 1 


- misl)Ff2 




5eV2 
2 


-P181)T^ 



''618,12618,/ 

''618,e618,4 " 



- Pi8l)Tj 
Pi8l)Ti 









^615 36,, ,={ex + Se{^-l)-n,,l)Zi 

' 15, a 

*^5,/-l5,.' = + ^'^"^ ""^'^^^a' 
tCye2, = iex + Se{4-l)-n^,l)ZI 

tl , =(6,+5e(^-l)-ni50^:; 

ir^.d 15, c 

tl^^^^,,^^,={ey-5,-m,,l)Zl, 

15, c 15, e 


te ,e,5 3 = i^y-^e-mi5l)F^ 

15, a ' 

,6 , =(e. + 3^-pi5/)T;, 

1 5 , a 1 5 , c 

*fi5,36^^_^, = (e?/ - ^6 - misO^s'^ 
e,_ , = (ej, - ^6 - "^isO^rf"; 

15. d lo,/ " 


^6^, ,.615,3 = (ex + Se{^ - 1) - ni5Z)F3 
,6,, , = {e^ + Sei^-l)-ni5l)F 



9) Vu = {ex + Se{^--^ + l),ey - 3|, e,) 





{v:,r 


{Vu)y 


— 

614,369^; 


(|e2,-3^| -muO^f 


fX 

6i4,jei4j 


= (6.- 


PiaI)T^ 


'-614,13614,; 


= {ez -PiaI)TI' 


ei4,feei4,j 




4.x 


= (e.- 


PlAl)Tf 


'■6l4,(.6i4,j 


= {ez-pul)Tj 


ei4,fcei4_; 


= (|ey-3^|-mi4Z)^f 


4X 

6l4,j6l4,a 


= (6.- 


Pul)Tf 






614,^614,13 


= (k,-3^|-mi4/)Zf3 


4X 

6l4,l36i4_fc 




-pi^m 






614,^614,13 


= (|ey-3^|-mi4/)Zf3 


614, 13614, /t 


= (k. 


-3||-mi4Z)Ff3 






ei4,fe6l4,a 


= {ex + 5e{^ + \)- nul)Zl 


6i4,fc6i4,;i 




-pul)T- 







Wc now need to define tlie off-diagonal elements of the matrix given by 



-1 



The terms t 



6i 



+ ^6, + represent the total number of stacks intersected by the edge Cj in the z-, x- and y-direction 
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respectively. Therefore we obtain 

= (5(sin(0e.) cos(6'eJ + sin(^ei) sin(6'ej + sin(^ei)) 

Throughout we wiU use the foUowing short-hand notation 



(5e(l + ^) for all edges such that 7^ 

(5(\/2 + 1) for those edges such that = 

(A.8) 



xvA-rii {\xv,\-nil) 



^(^/^/2 + l) ' 4(1 + ^ 



^, _ {\zVi\-Pil) ^ _ {\zVi\-Pil) 

Here, e^^j denotes the edge starting at vertex Vi and ending at vertex Vj 

We will now give the values of the upper-diagonal entries of the matrix ^/A ^ . 
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A 
A 

Va 
Va' 
Va 
Va 
Va 
Va 

A 
A 

Va 

A 
A 
A 

Va 
Va 
Va 
Va 

A 
A 

Va 
Va 
Va 
Va 
Va 
Va 

A 
A 
A 

Va 
Va 
Va 
Va 

A 



-1 



-1 



-1 



A 



h,i — 4 




, „ a' 

\lo,7o -A) 
~ 2 


{VA 


\2o,5o _ -Cq\/2 

/ K O A 

J 0,1 4 


{Va 


\02,32 -^2 

''9,2 


{Va 


\8o,7o — C'o 
/6,3 ~ 2 ^ 


{Va 


\6o,7o -^0 
''7,5 - 72 


{Va 


,,102,112 -C2\/2 


{Va 


ill, 7 ~ A 


\32,122 —B^ 


{Va 


J 10,9 — 72 


N,23,193 -A3 


{^ A 


''14,9 — ~2~ 


Nl22,ll2 -^2 

ill, 10 — 72 


{VA 


\ 153, 173 -Ca 


{Va 


n5,12 — 2 


\3l5,Ci5 -Cj^5\/2 

^39,12 — 4 


{Va 


\ 163,193 -B3 


{Va 


Jl4,13 - 72 




a73,133 _ -A3 


{Va 


^6,15 ~ 2 


\3i3,9i3 _ -C13V2 
/29,16 4 


{Va 


\a4,134 ""^4 

/ 1 n 1 — 
/iy,i/ 2 


{Va 


\184,134 -S4 
h9,18 - 72 


{Va 


N|4i8,ei8 -Ci8\/2 


{Va 


J32,18 — 4 


\4i8,ai8 —^4 


{Va 


J 36, 18 - 72 


X4l3,9l3 -b\^ 


{Va 


^29,19 — 2 


\C4,fe4 -A4 
J21,20 - 72 


{Va 


\65,P5 -As 
J24,23 — 2 


{Va 


^65,l35 -C5V2 
/'27,23 — 4 


{Va 


loo OA 


{Va 


\ 1^6,136 — 

^27,26 - 72 


{Va 


^5l3,9l3 -^13 

J29,27 - 72 


{Va 


\9i3,14i3 _ -Bi3 ( 2 ,-,\ 
^30,29 ~ 2 *>72~'~ '' 


{Va 


. 1 Q, . . f{' 


{Va 


M3,30 - 72 


^el8,/l8 _ -B18 (2 t -i\ 

)'\A'\2 — 9 \ /o ' ^ ) 
/ 0^,0^ A ^ ^2 


{Va 


\12l8,/l8 _ -C18V2 
/o4,oo 4 


{Va 


/ 

a/i8i9i8 _ -^18 
.'35,34 - ^ 

It t 


{Va 


/39,37 2 


(\fA 


N,(ij^5,C]^5 —^15 


{VA 


^39,38 ~ 2 


^' a' a' 


{Va 


Ml,40 - 72 


M6,42 — 2 


{Va 


\ai4,/H4 ""^14 





r-l 



I0180 _ 

3.1 ~ 

2(),3o 

3.2 — 

02,102 _ 

7,2 - 

02,112 _ 

11,2 - 

17o,6o _ 
5,4 - 

102,92 _ 

8.7 - 

92,32 _ 

9.8 ~ 

32,112 _ 

11,9 ~ 

23,173 _ 
15,9 - 

153,163 . 
13,12 

153,133 
16,12 

3l5,/l5 

40,12 " 

163,133 
16,13 

3l3,4l3 

19.16 " 

3l3,21l3 

31,16 

a4 ,C4 - 

20.17 ~ " 

184,64 _ 

21.18 ~ 

4l8,12l8 

33,18 

134, C4 

20.19 ~ 

4l3,14l3 

30,19 

£4,164 

22.20 ~ 

65,55 _ 

25.23 — " 

P6,?6 _ - 

25.24 — " 

96,135 _ 

27.25 — 

175,r-5 _ 

28.26 — 

5l3,14l3 _ 

30,27 

14l3,21i3 

31,30 

13l4,h4 _ 

44,30 

ei8,9i8 _ 
35,32 — 

12l8,9l8 

35,33 

/l8,lil8 _ 

36,34 — 

^15 'As 

40.37 ~ 

^'15 '/l5 _ 

40.38 ~ 

fcl4,'l4 

44,42 — 



-C0V2 



V2 
2 
2 

-Coa/2 



2 '^^^72'' 



2 

-C2V2 



2 

-B3 
V2 

-A3 
V2 

V2 
-C3V2 



Ha 



V2 



_ -Bu 
2 

-C4V2 



2 yV2^ > 



-B4 
2 



2 

^V2 

_B^ _ CbV2 
2 4 

-A 
2 



4 

— o,. 



_ Ms. 
2 

V2 



4 

_ -Ci3\/2 
4 

. — Ci4\/2 
4 



2 



2 



2 



<li43h4 

4§,'J3 ~ 

^\^14,jl4 



V2 

-C'i4a/2 
4 

Cl4V^ 



Va 
Va 
Va 
Va 

'A 
A 
A 
A 
A 
A 

Va 
Va 
Va 
Va 
Va 
Va 

A 
A 
A 

Va 
Va 
Va 
Va 
Va 

A 
A 
A 

Va 
Va 
Va 
Va 

A 



Bl4 



A 



-1 



-1 



lo,6o _ 

5.1 ~ 

2o,4o _ 

4.2 ~ 

02,92 _ 

8.2 - 

8o,17o _ 

4.3 - 

17o,7o _ 

6.4 - 

102,122 

10,7 

92,122 _ 

10.8 ~ 

23.163 _ 

13.9 ~ 

23,133 _ 
16,9 - 

153,193 
14,12 " 

3l5,dl5 

38,12 - 

3l5,<ll5 

41,12 - 

193,173 _ 
15,14 - 

3l3,5l3 _ 

27.16 " 
04,184 

18.17 ~ 

a4 ,64 

21.17 ~ 

184,164 

22.18 " 

4l8,9l8 _ 

35.18 - 

134,164 

22.19 ' 

4l3,21l3 

31,19 

64.164 

22,21 — 

65,175 _ 

26.23 — 

P5,175 _ 

26.24 - 

9s ,^6 _ 

28.25 — 

135 ,r5 _ 
28,27 — 

5l3,21l3 

31,27 

13i4,A;i4 
42,30 

13l4,il4 

46,30 

ei8,i^l8 _ 
36,32 - 

12l8,rfl8 

36,33 

e d' 
^^15 '"15 _ 

38,37 " 
/ / 

^15'"l6 _ 

41,37 " 
/ / 

•^15 '"^IS _ 

41,39 " 

fcl4,il4 
45,42 - 



C0V2 
4 

-Bn 
2 



Bo 
V2 



-C,V2 



-B3 
2 

2 

-B3 
2 



2 
2 

-B3 
2 



Al5 

V2 



2 

-B4 

2 V72 



-A4 
2 

"^^18 
2 

-A 
4 

2 '^72^ '' 



-B4 
V2 

-Cr,V2 
4 

z£k 
2 

4 



14 



2 

-Bm , 
2 ' 

-^18^ 
4 

-b' 



2 

V2 



+ 1) 



V2 

2 '^72^^'' 
2 



ai4,Jl4 
45,43 

^\jl4,^14 



2 

4 



A. 3 Edge-metric components for the 8-valent graph 



In this Section we give the values for the terms 
rotated, and translated graphs respectively. 



for the 8-valent graph as computed for aligned, 



A. 3.1 Aligned 8-valent graph 

For the case of an 8-valent graph, the periodicity cell contains five vertices, therefore \fA is a 36 x 36 
matrix. Below we give the non-zero upper-diagonal elements of the matrix \fA . The lower indices refer 
to the matrix elements while, the upper indices, indicate which edges we are taking into consideration. 



A 
T 
T 
T 

Va' 
Va' 

^/A~ 
^/A~ 

Va~ 

A~ 

T 
T 

A~ 



-1 


n0i,13i o 
'10,1 — P 




^0i,12i 

ai,i — 


2P 


f%/I"')°^'"l - - 


-B 

P 


-1 


n0i,9i lo 
/13,1 ~ 2^ 


(a/I"' 


^0i,15i 

a4,i — 


P 


IV^ Jl5^1 — 


2P 


-1 


a4i,13i la 
110,9 ~ 2^ 




a4i,i2i 
ai,9 — 


P 


Ta/I"')^^^'^^! - 


2P 


-1 


a4i,9i o 
h3,9 ~ P 




a4i,i5i 

a4,9 ~ 


— B 
P 


rA/I~')^^i'^°i - 


— B 

P 


-1 


^13l,12l Ir, 

hi, 10 ~ 2P 




asi.iii 
a2,io — 


— B 

P 


[y^ )i4,io — 


— B 
P 


-1 


^13i,10i la 
/15,10 ~ 2P 


v_ V 


a2i,iii 

'12,11 ~ 


2P 


(\/A)7},-, = -3 




-1 


/14,11 1-' 


i^/A^^ 
V V 


vl0i.9i 
a6,ii ~ 


-B 

p 


(^/A~^)ll''^T^^ - 
723,11 


-B 

P 


-1 


\ll2,ai2 no 

'24,11 P 




\li2,ei2 

'28,11 


-28 


(VAr')]\'f,' = - 

V » /13,12 


-iB 
2^ 


-1 


\lll,15i a 
h4,12 — ~P 




v9i,15i 
h4,13 — 


-¥ 


\^ )l5,13 — 


-P 


-1 


a6,13 — P 




a5i,ioi 

h5,14 — 


-¥ 


(Vl-')tl;^, = - 


P 


-1 


\59,Cg a 
^20,19 — P 




^22,21 — 


-p 


(a/I-')™ = 


-2/3 


-1 


\m-12,ei2 cyn 
/28,23 — ^P 


(a/I-^ 


\Cii2,ei2 
^28,24 — 


-a 


{y^ J26,25 — 


-2/3 


-1 


\bi2,hi2 a 
/27,25 ~ P 


{Va' 


\^'12.fcl2 

^29,2,5 — 


-w 


/ /T~^a4i2,/ii2 _ 
IV ^ J27,26 ~ 


-2/3 


-1 


\4i2,A;i2 a 
'29,26 — P 


{^~^ 


Jtl2,fcl2 

'29,27 ~ 


-2(5 


VV^ J31,30 - 


-P 


-1 


aie,14e 1 
/32,30 ~ 2P 




\lle,Pe 

'33,30 ~ 


-2/3 


i^'')llfo = - 


-\P 


-1 


\lle,ge _ 1 fD 

/36,30 ~ 2P 


(a/I-' 


a3e,14e 

'32,31 ~ 


-/3 


i^'')ll:fi = - 


-\P 


-1 


\13e,ge 

/34,31 ~ P 


(a/I-' 


\13e,re 

^35,31 ~ ■ 




IV ^ ^34 32 — 


4/3 


-1 


^14e,re _ r,o 

/35,32 — ■^P 


(a/I-' 


\14e,3e 

^36,32 — ' 


-¥ 


IV ^ ^34 33 — 


P 


-1 


\'Pe,re la 
/35,33 ~ 2P 


(a/I-' 


^36,33 ~ 


-p 


IV ^ j35_34 — 


P 


(V 


'~A-^\qe,ge _ lo 
^ /36,34 — 2P 


(VI-' 


^36,35 ~ 


-p 







where /3 



nl) 



5^3 



A.3.2 Rotated 8-Valent graph 

In this Section we give the values of all the terms that appear in the matrix \/I ' for the rotated 

8-valent graph. As in the aligned case, the matrix ^fA is a 36 x 36 matrix labelled by the edges of the five 
vertices that comprise the periodicity cell. The following short-hand notation is used: Hhi,kj ~ — 
which is defined in a similar manner as that done for the 4, and 6-valent graphs. 
The elements of the upper-half of the matrix \/I are: 
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2^12i,0i 




-1 

12,1 " 




— — //II n 

2A^lll,0l 


V ^ > 


-1 

10,1 " 




2^^131,01 


V ^;i3,i 




2^''9l,0l 




-1 

15,1 " 




~ 2^10ii0i 




-1 

15,1 " 




— i // 1 r n 

2P'15i,0i 


V ^^2,11 




~ 2^11i>12i 




-1 

10,9 " 




~ 2^13lil4i 




-1 

11,9 " 




~ 2^14i,12i 


'^/^)l3,ll 




~ 2^9i>12i 




-1 
15,11 




~ 2^10i,12i 




-1 
27,11 




— — // /I 1 

2H'4i2,ll2 


V ^^24, 11 




2^''12)ll2 




-1 
25,11 




— i ///i 1 

2PDl2,il2 




-1 
28,11 




2^'^12:ll2 


'^/^)27,11 




2^''nl2,ll2 




-1 
23,11 




— — //™ 1 

2/^mi2,li2 




-1 
12,10 




~ 2^13i,lli 


'^^'^)9,12 " 




~ 2^14i,lli 


(\fA: 


-1 
13,12 




1 ,, 

2/^91,111 




-1 
14,12 




2^15i,lli 






— — // 1 /I 

2M14e,ge 




-1 

10,9 " 




~ 2/^14i,13i 




-1 
15,10 




~ 2^10i,13i 


'^/^) 14,10 




~ 2^15i,13i 




-1 

13,9 " 




~ 2/^9lil4i 




-1 

15,9 " 




~ 2 ^lOl il4i 


V ^;i4,9 




~ 2^15iil4i 


v ^ > 


~1 
15,13 




~ 2^10ii9l 




-1 

14,15 




— 77 1 r n 

2^''15l,9l 


V^)l7,13 




— — // /I 1 




-1 
16,13 




2^(19,19 




-1 

18,13 




2^09>l9 


V ^721,13 




— — //ft 1 


(\fA: 


-1 
19,13 




— /Apr 1 

2^09,19 




-1 
22,13 




— —77^ T 

2 39,19 


V ^715,14 




~ 2^15i,10i 


v ^ > 


-1 
17,16 




— //^ A 

2^(19,49 




-1 

18,17 




~ 2 1^9 


^^)21,17 




— tin A 

2F'D9,49 


V ^ > 


-1 
19,17 




— // r A 

2P59,49 




-1 
20,17 




— — 77 „ A 

2^''C9,49 


V ^ns.ie 




2 09, 09 




-1 
21,16 




~ 2^'39,a9 


v ^ > 


-1 
22,16 




— — //j ^ 


V ^^20,16 




2H'Cg,ag 




-1 
19,18 




2^^59,09 




-1 
22,18 




2H'dg,bg 


V ^.'20,18 




— — // „ u 

2^''C9,6g 


(\fA: 

V ^ > 


-1 
21,19 




— — // cr a 

2^59,69 




-1 
22,21 




2 39,69 


V ^/'21,22 




2A^C9,fa9 


V ' > 


-1 
22,19 




2Prf9,59 




-1 
20,19 




~ 2^^9,59 


V ^^22,20 




2A*C9,(i9 


(\fA: 


-1 
32,31 




~ 2^13eil4e 




-1 
30,28 




~ 2^12e,lle 


V ^^33,30 




2^''Pe)lle 




-1 
34,30 




— — //« 11 

2A^ge,lle 




-1 
36,30 




— — // ^ 11 

2^36, lie 


V ^732,30 




~ 2^14e,lle 




-1 
31,30 




~ 2/^13e,lle 


(a/I^ 


-1 
26,24 




— 11 ^ A 

2A''ai2,4i2 






— i / / 7i /I 

2F'Ol2,4l2 




-1 
28,26 




1,, 

2^*612, 4l2 




-1 
27,26 




— I LU A 

2r'A!,12,4i2 


V ^729,26 




2Mfcl2,4i2 




-1 
25,24 




— // L 

2Pbi2,ai2 




-1 
28,24 




2P'ei2,ai2 


V ^^29, 33 




2^''ml2,al2 


V V 


-1 
29,24 




— — / / 7, 

2Pfci2,ai2 


(a/I^ 


-1 
27,25 




— — til. u 
2^/112,012 


V ^^25,23 




1 ,, 

2Mmi2,bi2 




-1 
29,25 




2Pbi2,«;i2 




-1 
28,27 




— — Ill, 

2P'All2,ei2 


V ^^28,23 




2Pmi2,ei2 




-1 
29,28 




— / / 7, 

2Pfci2,ei2 


(a/I^ 


-1 
35,28 




1 ,, 

2A^r-e,12e 


V ^733,28 




~ 2^P<!,12e 




-1 
36,23 




— — // ^ 10 

2P3e,12e 




-1 
32,28 




~ 2^14e,12e 






— ^Ml3„ 12„ 




-1 

07 OQ 
Z / ,ZO 




2r"tl2,'tl2 




-1 




-^/^fcl2 hl2 
2r^7^12i'tl2 


V^J29,23 




— 2^'=12:"^12 




-1 

35,33 




~ 2^Pe,re 


(^^; 


-1 
35,34 






V ^J36,35 




~ 2^9e,re 




-1 
35,32 




— 2^14e,re 


(^^; 


-1 
35,31 






V ^J34,33 




~ 2^96,^6 




-1 
36,33 




~ 2^9e,Pe 




-1 
33,31 




~ 2^13e,Pe 


V ^j36,33 




~ 2^9e,qe 




-1 
34,32 




— 2/^14e,ge 




-1 
34,31 




~ 2^13e,9e 



A. 3. 3 Translated 8-valent graph 

In this Section we will first give the values of all the terms te^e. = t^.^_. + t^^ej + 1^.^. as computed for all the 
five vertices comprising the periodicity cell for the 8-valent graph. We then give the values for the terms 

' J G ■ J I — 1 I — 1 

. that appear in the matrix \J A .As for the aligned case, the matrix \J A is a 36 x 36 matrix 

labelled by the edges of the five vertices that comprise the periodicity cell. 

The values for te^e are listed below. The (Gauss brackets) terms rij = [-|-], w-^"!"] and = [-|-] are 
the number of stacks that each edge intersects in the x-, y- and z-direction respectively (see Section [3. 1|) 

1) ^0" = (ea:,ej/,e^) 
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eo. 560,6 


= = 

60,560.7 




60,560,7 


= 

60,560,4 




''60,560,6 


''60,560,4 




60,860,6 


— — 

60,560,8 




^60,560,3 


~ ^60,460,7 




''60,460,6 


~ ^60,560,1 




60,360,1 


— 

60,760,6 




fX 

60,760,3 


— 

60,660,8 




fy 

''60,660,1 


— ''60,860,7 ~ 




60,460,3 


— — 
60,460,2 




60.660.1 


— 

60,660.2 




'■60,860,3 


-fy - 

— ''60,860,2 ~ 




60,260.1 


— f^ — 

60,160.4 


{ey - mol) 


60,860,1 


— — 
60,860,2 


{ez - Pol) 


fy 

'■60,760,3 


-fy - 

''60.760,2 


(e^ - Pol) 


60,360,2 


— /-^ — 
60,760,8 


(e^ - no/) 


60,160,2 


_ _ 
60,460,3 


(cy - mol) 


''60,260,3 


-fy - 

''60,460,1 


{^x - nol) 









f^ 


— f^ — 






fX 

''61.961.0 


61.1461.13 




-fy - 


61,1161.13 


''61,961.0 








"-61,961,10 — ''61.961,12 — 


V- 

61.1561.10 


— — 

61,961,15 






61.961.11 


— f^' 

61.1261,0 




''61,1261,10 — ''61,961,14 








f^ 


— 






fX 

61,061,11 


— 

61,1061,15 




"-61, 1261, 14 — '-61,1561,0 — 


61,1161,14 


61,061,10 








f^ 


— f^ — 






fX 


— f^ 




''61,1561,11 ''61,1361,11 


61,1261,11 


61,1261,13 






61,1261,11 


61,1061,13 




61,1361,14 


— — f 

61,061,15 V 




— mi^) 


61,1561,14 


_ j-a; _ 

61,1561,13 


(ley - ^1 - mil) 


- - 

''61,061,11 ''61,061,13 


61,961,10 


-t^ — ( 

61,1461,12 V 


~ 73I 


- ■«,!/) 


^61,961,12 


— f^ — ( 
61,1061,14 V 


1^^ ~ ^1 


= (ka.-^|-mO 












''61,1561,13 ''61,1061,14 

= (e. + ^ - mil) 



3)< = (e,,e,-2^,,^ + 5^) 





(^9> 




69.469,1 


— fZ 

■ 






fX 

69,069,5 


— f^ 

69.069.1 




fy 

'■69,c69.£j 


~ ''69,069,6 — 


f^ 
69, 669, d 


— f^ — 

69,c69,6 






^69,069,4 


— 

69,669,5 




fy 

'■69,669,d 


~ ^69,c69,a 


69,469, a 


— f^ 

69, 569, d 






fX 

69,569,4 


— f ^ 

69,^69,6 




fy 

''69,669,0 


''69,669,5 


f^ 

69,669,4 


— f ^ — 

69,669,1 






69.569,4 


= t'' 

e9,d69,i 




'■69,669,4 


- fy - 

— ''69,669,1 — 


69, 169, a 


— fZ _ 

69,569,6 






69, 669, a 


_ J-X _ 

69,669,1 


(|e,-2^|-m9Z) 


fy 

"-69, 569,4 


= ^69,569,1 = i^Z - P9I) 


69,c69,d 


— f^ — 

69,a69,6 




- ng/) 


69,069,6 


— fX _ 
eg,d69,a 


i^z - P9I) 


fy 

''69,d69,a 


— ^69,169,4 ~ {.^X 



4)^i'; = (e.,e,-^,6, + 2^) 









f^ 

612,6612,0 


— f^ — 

612,4612,1 




612,6612,4 


— f^ 

6l2,rn,ei2,e 






y?/ _ yJ/ _ 

''612.6612.0 ~ '■612,e6l2,;j 


f^ 

612, 1612, a 


_ +2 _ 

612,6612,1 




612,6612,(1 


— f^ 

ei2,fc612,4 






y?/ _ y?/ _ 

''ei2.fcei2.o '■612, 6612, m 


ei2,hei2,m 


— f^ 

612,4612, a 




y-a: 

"612,4612,;, 


— 

612,aei2,l 






y2/ _ y?/ _ 

''612, 0612, m ''612,1612,4 


f^ 

ei2,hei2,e 


— f^ — 
ei2,fcei2,e 




6l2,fc6i2,;i 


— P 

ei2,aei2,e 






y?/ _ yy _ 

''612,1612,^ ''612,1612,6 


f^ 

612,e612,m 


_ y-2 _ 

612,m6l2,fc 


{ex - ni2l) 


612, 1612, m. 


— f^ 

612,1612,6 




1 - mi2/) 


y!/ _ y!/ _ 

''6l2,fc6l2,m "-612,4612, 6 


6i2,fc6i2,h 


-t^ - ( 

612,6612,4 V 


kj/- 2^1 -mi 2/) 


ya; 

612, 6612, A: 


— f ^ 

612,a612,m, 


= (^^ + 4 


-Pi 2/) 


= (Ca; - '^l20 












y2/ _ yJ/ _ 

''612,6612,*: '■612,4612,(1 

(e, + 2^-;>i2/) 



C\ T/" _ SeV2 , _ SeV2 ^ I q W2\ 

^^6 — l^^a; 2 ' ^2/ 2 ' ^2 "-^ 2 
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iv"^y 

\^ e ) 




— 

ee,l2ee,ll 








ee,12ee,ll 


^6e,r6e,g 






fV -fV _ 

Ge.l2Ge.l4 Ce.l2 6e,g 
















— +2 

ee,12ee,13 








4.x 

ee,12ee,p 


= 

66,366,11 






- f - 

6e.g6e,14 G.e 12^e r — 


ee,13ee,14 












*e e = 












= 

Cg 14Ce 13 






fV ' -fV ' ' - 

^G-r- n€.e r '^€e I'^Cc 11 


llCe 14 














' ' = 


' ' = 
















''ee,1366,p ~ ''6e,q6e,p — 












66, 1466, q 






= 


— 

ee,liee,13 


(l%- 


5 1 
f 


- rrieO 




~ ^6e, 1366,9 ~ 




- 


''66, 1166, p — ''6e,ll6e,q " 


fZ 

Ee, 1266, 14 


= = 






- Hel) 


''66,1266,9 


= t^ = 
Ce,14Ce,r 




-Pel) 


















''66,1366,9 "-66,1466,7- 






































= (e. + 3^-PeO 



The values for the upper-half of the matrix \YA are given in the following table. As for the 4- and 6- 
valent graphs, the lower indices represent the matrix entries while the upper indices indicate the edges under 
consideration. 
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\5(),6o 
12,1 — 




1 

iVA-' 


\5o,7o 
'3,1 ~ 


1^0 


1 

{VA-' 


\5o,8o 
'4,1 ~ 


-\Bo 




\5o,4o _ 
'5,1 ~ 


Co 


, 1 

{VA-' 


\5o,3o 

^6,1 ~ 


\Co 


, 1 

{Va-' 


\5o,lo _ 
'8,1 ~ 


-\Co 




\6o,7o 
l3,2 — 


I Bo 


, 1 


\6o,8o 

^4,2 — 


'2D0 


, 1 

{VA-' 


\6o,4o 
'5,2 - 


-\Co 


[Vat' 


\6o,2o 
17,2 — 


^Co 


iVA-' 


\6o,lo 
^8,2 — 


2Co 


{VA-' 


\7o,8o 
'4,3 ~ 


-\Eo 




JoAo 
15,3 — 


^Co 


iVA-' 


\7o,3o 
^6,3 ~ 


Co 


{Va-' 


\7(),2o 
'7,3 ~ 


-Co 




\8o,3o 
'6,4 ~ 


^Co 


iVA-' 


\8o,2o _ 
^7,4 — 


Co 


{VA-' 


\8o,lo _ 
'8,4 - 


-\Co 




\4o,3o 
^6,5 ~ 


Bo 


{Va~^ 


n4o,2o 
^7,5 — 


^Bo 


{Va'' 


\4o,lo 
'8,5 ~ 


-\Fo 


[Vat' 


\3o,2o 
>7,6 - 


Ao 


{Va.-^ 


\3o,lo 
^8,6 ~ 


^Bq 


{Va.-' 


\2o,lo 
'8,7 - 


-Bo 


[Va 


v0i,9i 
^9,8 ~ 


Bi 


1 

{Va~^ 


\0i,10i 
/10,8 ~ ■ 


-\Bi 


{Va-' 


\0i,15i 
'11,8 ~ 


---Fx 


[Vat' 


v0i,12i 
a3,8 ~ " 




{Va~^ 


^Ol,lll 
/14,8 ~ ■ 


-\F, 


{Va.-' 


n0i,13i 
'15,8 ~ 


--\A, 


[Va 


v9i,10i 
ao,9 ~ " 


-Ai 


1 

{Va~^ 


\9i,15i 
/11,9 ~ ■ 


-2-Bi 


{Va.-' 


\9i,14i 
'12,9 ~ 


--'2M 


[y/A-' 


\9i,12i 
a3,9 ~ " 


-A, 




\9i,lli 
a4,9 ~ ■ 




, 1 

{VA-' 


\l9,cf9 

'l6,9 ~ 


-\Bg 




\l9,59 

a8,9 — 


^Cg 


1 


\ 19,69 

a9,9 — 


\Dg 


1 

{Va-' 


a9,49 
'19,9 — 


-\.Fg 


[^^A 


\l9,69 

^20,9 — 




1 

{^^A 


\l9,a9 
^21,9 — 


-Bg 


{Va-' 


\10i,15i 

'11,10 


= -Bi 




aoi,i4i 
a2,io ~ 


-\Ei 




\10i,12i 

a3,io ~ 




{Va-' 


\10i,13i 
'15,10 


= -^Bi 




a5i,i4i 

a2,ii — 


-\Bi 


1 

i^fA-' 


a5i,lli 
^4,11 — 


-1^1 


1 

{Va.-' 


sl5i,13i 
'15,11 


= --2^1 


, 1 

[VA-' 


\14i,12i 

a3,i2 — 


-Ai 


, 1 

{VA-' 


\14i,lli 
a4,12 — 


-\Bi 


{Va-' 


\14i,13i 
'15,12 


= -Bi 


[Va~' 


a2i,iii 

a4,i3 ~ 


-Bi 


1 

iVA-' 


a2i,13i 
a5,13 ~ 


-\Bi 


1 

{Va-' 


\li2,ai2 
'23,13 " 


= -\Fi2 


, 1 

[VA 


\ll2,6l2 
'24,13 ~ 


-1^2 


, 1 

{Va~' 


\ll2,4l2 

^25,13 ~ 


-A12 


, 1 

{VA-' 


\li2,mi2 
'26,13 


= ~\Bi2 


[Va ' 


\ll2,'ll2 

^28,13 — 


-\M2 


, 1 

{Va-' 


\li2,ei2 
^29,13 ~ 


-\Fi2 


, 1 

{Va-' 


\lli,13i 
'15, 14 


= -\F^ 




^d<.),cg 

h5,16 — ~ 


\Bg 


{Va.-' 


\<i9,59 

a8,16 ~ 


-\Bg 


{Va.-' 


\rf9,69 

'19, 16 ~ 


-Bg 


1 

[^/A~^ 


\d9,bg 
^21,16 ~ 




{Vat' 


>(i9,a9 
^22, 16 ~ 


-\Eg 


{Vat' 


\C9,59 _ 

'l8,17 ~ 


-Bg 


[^/A~^ 


\c<j,6<j 

a9,i7 ~ ~ 


\Bg 


, 1 

{Vat' 


\C9,49 

^20,17 ~ 


\Bg 


, 1 

{Va ' 


\C9.69 
'21,17 ~ 


-Cg 


[^^Ar^ 


\C9,09 

^22,17 ~ 


\c% 


{Vat' 


\59,69 _ 

/19,18 ~ 


-hDg 


{Vat' 


\59,49 _ 

'20,18 ~ 




[^^A 


\59,69 

^21,18 ~ 

\ 49,69 
^21,20 — 


\Bg 


{Vat' 


\69,49 

/20,19 ~ 

\49,a9 
^22,20 — 


-\Cg 


{Vat' 


\69,«9 
'22,19 ~ 
\69,09 
'22,21 — 


-\Bg 


[^^A-^ 


Bg 


{VT' 


-\Bg 


{VT' 


-\Eg 


, 1 

{Va~' 


vai2,6i2 
^24,23 — 


-A12 


{VT' 


,ai2,4i2 
^25,23 — 


-hA,2 


{Va.-' 


\ai2,mi2 
'26,23 


= -\Ei2 




,ai2,fci2 
^27,23 — 


-\M2 


{VT' 


\ai2,ei2 
^29,23 — 


-5-B12 


{VA-' 


\6l2,14l2 
'25,24 


= —B12 


[^-' 


\6i2,mi2 
^26,24 — 


-\M2 


1 

{VT' 


\6l2,fcl2 
^27,24 — 


— C12 


{VA-' 


\6l2,/ll2 
'28,24 ■ 


= -\Bi2 




a4i2,A:i2 
^27,25 


-\Bvi 


1 

{VA-' 


\14l2,/ll2 
^28,25 


= -1^12 


1 

{Va-' 


\14i2,ei2 
'29,25 


= -1^12 


[VA-' 


vmi2,A;i2 
^27,26 ~ 


-Avi 


1 

{VA-' 


>mi2,/ii2 
^28,26 ~ 


-1^12 


1 

{VA-' 


\mi2,ei2 
'29,26 


= -\Fi2 




,fcl2,/ll2 

^28,27 — 


—B12 


{VT' 


\fci2,ei2 
^29,27 ~ 


-1^12 


{VT' 


\fti2,ei2 
'29,28 " 


= -A12 


[Va 


^12e,lle 

^30,29 ~ 


-Be 


, 1 

{VA-' 


,12e,14e 

'31,29 ~ 


-Ae 


{VA-' 


\12e,13e 

'32,29 


= -\Be 




\12e,ge 

^33,29 ~ " 




, 1 

{VT' 


\12e,Pe 
'34,29 ~ ■ 


-^Be 


, 1 

{Va-' 


\12e,re 

'36,29 ~ 


2^e 


[Va ' 


,lle,14e 

^31,30 — 


-\Be 


{Va-' 


\lle,13e 
'32,30 ~ 




, 1 

{Va-' 


\lle,5e 
'33,30 ~ 


-- -\Be 


[VaT^ 


\lle,Pe 

^34,30 ~ " 




{Va ' 


\lle,'7e 
'35,30 ~ ■ 


2^e 


{Va ' 


\14e,13e 
'32,31 


= -Be 


, — 1 
[Va 


\14e,ffe 

^33,31 ~ " 


-1/4 


1 — 1 
{Va 


1/1 ^ 
'35,31 ~ " 


-\Be 


1 — 1 
{Va 


1/1 rr. 

'36,31 ~ 


---Ce 


[Vat' 


\13e,Pe 

^34,32 — " 




{Va~^ 


\13e,ge 
'35,32 — " 


-kDe 


{Va-^ 


\13e,re 

'36,32 - 


-- -\Be 


[Vat^ 


\9e,Pe _ 

1 34,33 — 


Be 


{Vat' 


\9e,«e _ 
'35,33 ~ 


-\Be 


{Va-' 


\fle,r-e _ 

'36,33 ~ 


-Ae 


[Va~' 


\Pe,qe _ 

^35,34 — 


If 


{Va~' 


\Pe,r-e _ 

'36,34 — 


-\Be 


{Va-' 


\9e,re _ 

'36,35 ~ 


-Be 



68 



where = (xv^-riil)^, Bi = {yy^-mil)-^, d = [zy^-pil)-^, Di = Bi + d, Ei = Ai+d, Fi = Bi+Ai, 
and the quantities xy^, yv^, zy^ are the x-, y-, z-coordinates of the vertex Vi. In ah the above formulae the 
term j^j^ comes from the fact that te^ = t%. + te^ + 1^. = ^eV^ for all edges Cj. 
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